N73-29885 



(L1SA-CR-1243N7) PREDICTION OF SPACE 
SHUTTLE FLUCTUATING PBESSOBE EN7IBONHENTS, 

INCLUDING SOCKET PLUME EFFECTS Final 
Report (Byie Labs., Inc.) 156 p HC Unclas 

$10.00 CSCL 22B G3/31 15228 












WYLE LABORATORIES - RESEARCH STAFF 
REPORT WR 73-6 


PREDICTION OF SPACE SHUTTLE 
FLUCTUATING PRESSURE ENVIRONMENTS, 
INCLUDING ROCKET PLUME EFFECTS 


By 

Kenneth J. Plctkin 
J. t . Rober son 


June 5973 


FINAL REPORT 
of 

Work Performed Under Cortroct NASO-26919 


ABSTRACT 


Preliminary estimates of Space Shuttle fluctuating pressure environment? have been made based 
on prediction techniques developed by Wyle Laboratories. Particular emphasis has been given 
to the transonic speed regime during launch of a parallel-burn Space Shuttle configuration. 

A baseline configuration consisting of a lightweight orblter and monolithic SR3, together with 
a typical flight trajectory, have been used as models for the predictions. Critical fluctuating 
pressure environments are predated at transonic Mach numbers with the following maximum 
overall fluctuating pressure levels: 1) shoulder-induced flow separation (157 dB), 2) terminal 
shock-wave oscillation (164 dB), 3) protuberance flows (168 dB), 4) compression-corner 
induced separated flow (148-157 dB), 5) shock-induced interference flow (161 dB), 

6) shock-free interference flow (150 dB), and 7) attached turbulent boundary layer flow 
(137 dB). Comparisons between predicted environments and wind tunnei test resists, in general, 
showed good agreement. Predicted one-third octave bend spectra for the above environments 
were generally one of three types: 1) attcched turbulent boundary layer spectra (typically 
high frequencies), 2) homogeneous separated flow and shock-free interference flow spectra 
(typically intermediate frequencies), and 3) shock -oscillation and shock-induced interference 
flow spectra (typically low frequencies). Predictions of plume induced separated flow environ- 
ments were made. Only the SRB plumes are important, with fluctuating levels comparable to 
compression-corner induced separated flow and shock oscillation. 
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1 .0 INTRODUCTION 


During launch and flight through the atmosphere, the external surfaces of the Space Shuttle 
vehicles will be exposed to high Intensity acoustic and fluctuating pressure environments. 

These environments will be critical to the design and success of the complete Space Shuttle 
system. Definition of the environments Is necessary, therefore, in order to insure structural 
integrity, reliability and the economic operational requirements of the Space Shuttle. 

In order to arrive at an optimized Space Shuttle system, accurate definition of the critical 
acoustic environments which will be encountered during Space Shuttle flight is required at 
the earliest possible time during the development stage. 

Wyle Laboratories is engaged in a research program under NASA-MSFC Contract NAS8-26919 
to develop reliable methods of predicting critical acoustic environment ■‘•'ends during Space 
Shuttle flight. The primary objective of this study is to o* Srm a critical analysis of the 
acousric and fluctuating pressure environments for the c lete flight envelope of rhe Space 
Shuttle. Particular emphasis is to be directed to the area of fluctuaring pressure environments 
as encountered during flight through the atmosphere. Key features of this investigation have 
been the bringing together of data f om disparate sources, combining the results to show general 
trends for the various unsteady phenomena, the development of prediction methods based on the 
combined results, and the application of the prediction methods in the final assessment of 
critical environments for the Space Shuttle vehicle. 

In early phases of this program, effort was devoted to a study of basic aerodynamic and acoustic 
environments for the purpose of: 1) formulating prediction methods for the basic fluctuating 
pressure environments which are encountered by aerospace vehicles; and 2) definition of 
the fluctuating pressure environments which are unique to the Space Shuttle vehicles. During 
this phase of the study, emphasis was given to the transonic speed regime during launch. 

During the past year, work was continued on the development of prediction techniques. 
Prediction schemes were refined and incorporated into computer programs. 

During the past year, flow field analysis and fluctuating environment prediction schemes were 
extended to include effects of rocket exhaust plume-induced separation. An extensive 
investigation of plume-induced separation was undertaken for the purpose of: 1) developing 
a scheme to predict the extent of plume-induced separation regions, and 2) evaluating the 
fluctuating environment associated with such regions. The results of this study were the 
subject of a recent report, and are incorporated into the predictions of the present report. 

This report presents the latest predictions for a Space Shuttle launch configuration. Selection 
of the vehicle booster and orbiter design and specification of trajectory parameters is presently 
in progress. Although these vital criteria have not been finalized, it is necessary to proceed 
with engineering studies to arrive at preliminary estimates of both steady and unsteady aero- 
dynamic loads which will be encountered during atmospheric flight bv the Space Shuttle 
vehicle. This report provides a cursory analysis of the fluctuating pressure environments which 
may be anticipated during launch, and predicted fluctuating pressure levels and spectra for 
critical flight conditions. These predictions are based on a typical Space Shuttle configuration 
and launch trajectory, and therefore, should be considered as preliminary estimates. 
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SecHon 2.0 contains a brief description of various types of fluctuating pressure environments 
which may be experienced by the Space Shuttle vehicle. Section 3.0 presents a description 
of the prediction methods developed by Wyle Laboratories for various fluctuating pressure 
environments together with comparisons of the predictions with experimental data. A 
description of anticipated unsteady flow fields and attendant fluctuating pressure environments 
which will be experienced by the mated Space Shuttle launch configuration is presented in 
Section 4.0. Section 5.0 presents an assessment of critical Space Shuttle flight conditions 
and anticipated overall pressure levels and spectra. Particular emphasis is given to the 
transonic Mach number range and the condition of maximum dynamic pressure. Finally, 
Section 6.0 summarizes the more significant conclusions reached during this stud/. 
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2.0 FLUCTUATING PRESSURE ENVIRONMENTS 


When a vehicle moves through air there are two basic means by which it can produce noise: 

1) by its propulsion mechanisms (motor-jet, rocket, etc.,) and 2) by i*s interaction with its 
surroundings. At low speeds, for example, during and immediately after lift-off, the first of 
these is by far the dominant one, while near or above rhe speed of sound the second becomes 
important. 

During any flight cycle for an aerospace vehicle, there are four important phases of the flight 
which should be investigated in order to assess the environmental trends due to acoustic and 
fluctuating pressures. These are listed in the chronological order in which they cc^ur. 

• Lift-off phase during which acoustic excitation results 
from the exhaust noise. 

• Launch Might to orbit phase, during which rocket exhaust noise 
diminishes and aerodynamic fluctuating pressures (pseudo-sound) 
start to dominate. From an aerodynamic noise viewpoint, this 
phase becomes most critical at transonic Mach numbers 

(0.60 < M <1.6). 

• Re-entry phase during which only aerodynamic fluctuating 
pressures are present. 

• Flyback phase during which Mie noise from flyback engines 
dominate. 

This report is devoted to the specification of surface fluctuating pressures resulting from unsteady 
aerodynamic flows which occur during the launch phase of flight since this is considered to be 
the primary area of study for this contractual effort. Inflight surface fluctuating pressures are 
distinctly different from acoustic disturbances originating from rocket exhaust flows and engine 
noise. Inflight disturbances, as considered herein, arise from several modes of disturbances — 
the principle source being the passage of a turbulence environment over the external surface. 
One other important source to be discussed is shock -wave oscillation which is characterized by 
both turbulence (in close proximity to the foot of the shock wave) and pseudo-static disturbances 
resulting from the modulation of the pressure gradient through the shock wave. Thus, inflight 
fluctuating pressure phenomena are near-field mechanisms acting on the surface of a vehicle 
with the distinction that the disturbances are generally convec ted at some fraction of the locai 
mean flow velocity. On the other hand, rocket noise and engine noise are acoustic mode 
disturbances which generally originate away from the surface. Furthermore, acoustic mode 
disturbances consist of sound waves which propagate at the local speed of sound with a direction 
independent of the local velocity. 

Two specific areas crs considered in the present study: 1 ) "basic" fluctuating pressure environ- 
ments which will occur on virtually all aerospace vehicles during some phase of flight, and 

2) "special" fluctuating pressure environments which are unique to the Space Shuttle 
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configuration and mission parameters. Under the first category fall such environments as 
attached turbulent boundary layers / separated flows and certain shock-oscillation environments. 
Other environments which are more dependent on configuration and mission parameters are 
certain types of shock-wave impingement, protuberance flows, wake and base flows, rocket 
plume interference, and rocket exhaust impingement. Distinction is made between the two 
areas since the "basic" fluctuating pressure environments are dictated by the aerodynamic 
flow field of the overall vehicle; whereas, the "special" fluctuating pressure environments 
are relatively localized flow fields which are essentially superimposed upon the basic flow 
field. 

There are an infinite number of possible vehicle configurations and any discussion of their 
fluctuating pressure environments must be general. Even simple vehicle shapes, such as 
cone-cylinder bodies, produce complex and highly nonhomogeneous unsteady flow fields — 
particularly at subsonic and transonic Mach numbers. Unsteady flow phenomena are of 
particular importance at transom* speeds, since in this range, fluctuating pressures reach 
maximum levels due to their proportionality to dynamic pressure (q). However, peak fluctuating 
pressures do not necessarily occu at maximum dynamic pressure for certain regions of a vehicle 
due to the nonhomogeneous nature of the flow field. For example, regions of the vehicle 
exposed to separated flow and the impingement of osci Mating-shock waves will experience 
fluctuating pressures at least an order of magnitude greater than regions exposed to attached 
flow. Thus, if separated flow and osciilating-shock waves are present, say at Mach numbers 
other than the range of maximum dynamic pressure, then peak fluctuating pressures will also be 
encountered at conditions other than at maximum q. Thus, it is easily seen that vehicle con- 
figuration is very important in the specification of fluctuating pressure levels since the source 
phenomena are highly configuration dependent in addition to varying with Mach number and 
angle of attack. 

In light of the foregoing discussion, one general statement can be made in regard to aero- 
dynamic fluctuating pressures. Regions exposed to the :me unsteady phenomenon will 
experience fluctuating pressure levels which are proportional to dynamic pressure. Thus, it 
can be readily seen that a fundamental parameter in the specification of the surface excitation 
is free-stream dynamic pressure and its variation with Mach number. For a given configuration, 
Mach number and angle of artack define the phenomena, and dynamic pressure defines the 
fluctuating pressure levels associated with the phenomena. 

In order to assess the fluctuating pressure environments for a vehicle of arbitrary geometry, such 
as the Space Shuttle vehicle, it is convenient to develop prediction formula for the statistical 
properties of the fluctuating pressures for each basic type of unsteady flow field. From previous 
wind tunnel tests and flight data, the following unsteady flow fields have been identified as 
sources of fluctuating pressures: 

a) Attached t^“hulent boundary layers 

b) Separated flows 

1 ) Expansion Corners 

2) Compression Corners 
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c) Oscillating shock waves 

1) Transonic terminal shock 

2) Supersonic compression shock 

3) Supersonic impinging shock 

d) Wake and base flows 

e) Jet/rocket exhaust flow interaction with the vehicle 

aerodynamic flow field 

f) Protuberance flows 

g) Rocket exhaust impingement 

During past research programs, not all o* these phenomena have been studied in sufficient 
detail to enable prediction methods to be developed (for example, items d and g). However, 
the other flow fields have been examined sufficiently to provide at least a first order analysis 
of the attendant fluctuating pressure environments. 

The surface fluctuating pressures for eac.> of the above flow conditions may be expected to 
exhibit different statistical characteristics due to dissimilarities in the mechanisms at work. 

For example, attached flow pressure fluctuations result from the disturbances within turbulent 
boundary layers. Separated flow pressure fluctuations result from disturbances within the 
separated shear layer and instabilities associated with the separation and reattachment points. 
Pressure fluctuations for shock-wave oscillation result from the movement of the shock wave 
and the static pressure discontinuity associated with the shock wave. Generally, shock-wave 
impingement occurs in ♦■he presence of either attached or separated flow and added disturbances 
result due to the combined environments. Thus, in order to develop prediction techniques for 
the many unsteady flow fields which a vehicle mav encounter, it is necessary to examine each 
flow field independently, and to define, at least empirically, the statistical properties of the 
attendant fluctuating pressure environment. 

The statisticai characteristics of each fluctuating pressure environment that are important in 
the analysis of structural response may be classified under three parameters: 

• The overall fluctuating pressure level 

• The power spectra 

• The cross-power spectra (or narrow band cross-correlation) 

In the following section, a summary of the prediction formula is presented together with typical 
comparisons of the predictions with existing experimental data. 
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3.0 GENERAL PREDICTION FORMULAE 


Analyses of various unsteady flow environments have been performed under th'$ contract for 
the purpose of defining empirical formulae for estimating the statistical properties of the 
attendant fluctuating pressures. This work is based on the results contained in References 
1 through 25, and the findings are presented in References 26-29, A detailed summary of 
the basic findings and prediction formulae is presented in Appendix A, The surface 
fluctuating pressures for both "basic" unsteady flows and certain "special" unsteady flows 
unique to the Space Shuttle vehicles have been reviewed. The resulting empirical formulae 
have been generalized for the various environments : n order to achieve 1) simplicity in 
their mathematical formulation, 2) a relationship to the physical parameters of the flow, 
and 3) good agreement with the available experimental data. As noted in the previous 
section, complete definition of a fluctuating pressure environment requires that the overall 
level, power spectra, and cross-power spectra be defined. A brief summary of prediction 
methods for each of these statistical parameters is given in the following subsections for 
attached flow, separated flow, and shock-wave oscillation. 

3.1 Overall Fluctuating Pressure Level 


The correct method of presenting overall fluctuating pressure levels for surfaces beneath the 
convected turbulence in boundary layers is in terms of the root-mean-square fluctuating 

pressure level, \ . Free-stream dynamic pressure, , local dynamic pressure, q , 

I — x ~ i 

and wall shear stress, t , have been used to normalize \ P 2 so that meaningful data 

cj 

collapse can be realized throughout the Mach number range* The most generally accepted 
normalizing parameter is and, thus, it will be used in the current expressions. 


The effects of Mach number on the normalized rm$ intensities of the fluctuating pressures in 
attached flows are shown in Figure 1 , There is significant scatter in the data which may be 
attributed to several factors: 1) background noi c e and free-stream turbulence in the testing 
medium; 2) instrumentation quality and the precision of the experimental technique, and 
3) data acquisition and reduction techniques, etc. For the range of Mach numbers covered 
in the data of Figure 1 , the normalized rms value of the fluctuating pressure varies from 




00 


0.006 at subsonic Mach numbers to 0,002 at supersonic Mach numbers. 


Lowson, Reference 1 , proposed the following semi-empirical prediction formula which appears 
to agree with the general trend in the data: 



0. 006/(1 + 0.14 M = ) 


( 1 ) 


It i* important to note that this formulc has some theoretical basis and is not strictly an 
empirical approximation of measured results (see Reference 1 ). The use of this formula at high 
supersonic and hypersonic Mach numbers should be done so with the understanding that it has 
not been verified in this Mach number range and may lead to significant error. Predictions for 


6 



high supersonic and hypersonic Mach numbers have been developed under this program and 
are presented in Reference 27. However, in the Mach number range up to, say M =3.0, 
the above equation is in good agreement with experimental results. 

It should be noted ihat the results presented in Figure 1 , particularly the wind tunnel result's, 
were obtained for both homogeneous and stationary flows at free-stream conditions and in the 
absence of external pressure gradients. Consideration should be given to local conditions 
which deviate from free-stream conditions. This work is currently in progress. 

The variation of fluctuating pressure level, norma red by free-stream dynamic pressure, with 
local Mach number for various separated flow environments is presented in Figure 2. These 
data represent both expansion and compression corners. For expansion corners, the largest 
levels occurred at low Mach numbers and decreased as local Mach number increased. These 
data represent the region of plateau-static pressure and the tolerance brockets on the data 
represent the variations in fluctuating pressure level within the region of constant static 
pressu.e rather than scatter in the measurements. A good empirical approximation to the 
expansion induced separated flow disturbances is: 



This equation is similar in form with that previously proposed for attached turbulent boundary 
layers. 

Fluctuating pressure measurements for the region of plateau-static pressure upstream of com- 
pression corners exhibited a somewhat different trend with Mach number. In general, the 
compression corner data showed un increase in fluctuating pressure level with increasing free- 

stream Mach numbei in the ran < , 1 .0 < M < 2.0 — reaching a constant level at Mach 

oo 

numbers ab^ "2.0. Free-stream Mach number is used here because adequate data is not 
available for determining the local Mach numoer in the vicinity of the compression-induced 
separated flow region. Derivation of an empirical prediction fomnula for the fluctuating pressure 
level within compression-induced separated flows will require further study. 

Generally, shock-wave oscillation produces the most intense fluctuating pressure levels that 
are usually encountered by a vehicle. Typical shock waves encountered by vehicles are: 

• Terminal shock waves for regions of transonic flow 

• Displaced oblique shock waves as induced by the separated 
flow in compression corners at local supersonic speeds 
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• Reattachment shock waves in the vicinity of the reattachmerit point For 
separated flows generated by both compression and expansion corners 

• Impingement shock waves as caused by local bodies such as strap-on 
rockets 

» Oblique shock waves associated with separated flow induced by highly 

expanded rocket plume at high altitudes. 

All shock waves may be expected to produce generally similar fluctuating pressure environments 
since the pressure fluctuations arise from motion of the shock wave. The mechanism of this 
motion was investigated under this contract, and it was found that the oscillation ~rz uiiven 
by veiocity fluctuations in the incoming turbulent boundary layer. T he maximum displacement 
is governed by mean flow conditions. Except for alternating flow conditions, where oscillation 
distance is much larger than shock thickness and the power spectrum is therefore a Poisson 
distribution, spectra and overall levels for various shock environments are quite similar. A 
discussion of shock oscillation and results of analytical investigation are presented In 
Reference 28, and are summarized in Appendix B. The overall 'evel is a strong function of 
Mach number in the transonic regime; typical overall levels for terminal shock wave 
oscillation are presented in Figure 3. 

3.2 Power Spectra 

Power spectra represent the distributions of the mean -square ^actuating pressure with frequency. 
Power speclra for the various fluctuating pressure environments are found to scale on a Strouhal 
number basis; that is, the frequency is normalized by multiplying by a typical length and 
dividing by a typical velocity. Normalized spectra enable similar, stationary flows to be 
represented by c single spectrum regardless of the scale of the flow field or the free-stream 
velocity. 

Numerous studies have been conducted to determine the proper parameters to be used to 
nondimensionalize the spectra for various aero-acoustic environments. Unfortunately, the 
choice of parameters which best collapses the data appear *o be dependent on the nature of 
the flucti atina pressure environment. In general, free-stream velocity is used as the normaliz- 
ing velocity parameter, although a typical eddy convection velocity (itself a function of 
frequency) has been used occasionally . The local convection velocity appears to correspond 
more closely with the physical situation for fluctuating pressures due to turbulent eddies. 
Selection of a typical length is more difficult. Boundary layer thickness, 6 , displacement 

thickness, 6* , wall shear stress, t , and momentum thickness, 8, have all been used by 

u 

various investigators for attached flow. For separated flow and shock-wave oscillation, local 
boundary layer thickness and separation length hove also been used. 

Generally, prior knowledge of the overall fluctuating pressure level is required to predict the 
power spectra; however, prediction formula for the overall levels have been developed to a 
limited degree as discussed in Section 3.1 . 
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Prediction formulae for the power spectra for the various unsteady flows, with the exception of 
shock -•wave oscillation, ha^e the following form: 


0(f) U 




:3> 


where 

u 

= characteristic velocity 


L 

= characteristic length 


f 

0 

=■ characteristic frequency 


n, k 

= spectrum shape factor 


Some mathematical properties of Equation (3), in particular the relationship between n, k 
and overall level , are discussed in Appendix C. 

For shock-wave oscillation, the power spectrum is a combination of power spectra for inviscid 
shock-wave motion as well as contributions from the separated flow near the foot of the shock 
wave. The prediction formula has the following form: 



where k^ m 0.25, the subscripts SW and S denote shock wave and separated flow 

respectively, and the superscript ! denotes the absence of viscous effects. A complete 
discussion of the power spectra predictions for the various fluctuating pressure environments are 
oresented in Reference 26. Comparisons of predicted power spectra with experimental data are 
^ ented for attached flow, separated flow, and shock-wave oscillation in Figures 4, 5 and 6. 
Also, power spectra predictions have been developed for nonhomogeneous attached and 
separated flows (see Reference 26). Comparisons of these predictions with experimental data 
are shewn in Figure 7 together with attached flow, separated flow and shock -wove oscillation. 
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3*3 Cross- Power Spectra 


The final requirement in determining the characteristics of the fluctuating pressure environment 
of an unsteady flow field is to define the narrow band, space correlation function or co-power 
spectral density. This parameter is the key function needed to describe on impinging pressure 
field on a structure in order to calculate the induced mean-square response of the structure 
(see, for example. Reference 17 for the structural response computational technique). The 
spatial correlation properties of a fluctuating pressure field can be obtained only from a careful 
and detailed examination of the flow field at a large number of points. 

Measurements by several investigators have shown that the co-power spectral density of the 
turbulent boundary layer pressure fluctuations in the direction of the flow can be approximated 
by on exponentially damped cosine function, and the lateral co-spectrcl density can be 
approximated by an exponential function. 

The spatial coherence of the fluctuating pressure environments defined by the co-power 
spectra have been evaluated for the various unsteady flows. The cross -power spectra for 
attached turbulent boundary layers, two- end three-dimensional separated flows and protu- 
berance-induced wake flows appear to be fairly well defined at least in the longitudinal 
direction. Much uncertainty remai -is for transverse cross-power spectra in upstream separated 
fiows, and protuberance wake flows. Both longitudinal and transverse cross-power spectra for 
the regions beneath oscillating shock waves require additional study. 

In general, the normalized co-power spectra for the various fluctuating pressure environments 
can be defined as exponentially darrped sinusoids for the longitudinal axis as follows: 


cu,o = «*Jii (5> 


Along the transverse ax>$, the normalized co-spectra may be defined by an e* ponential decay 
as follows: 


C(n,f) = e" bun/U c (6) 

where u is the circulo* frequency, is the convection velocity, and £ and n are the 

longitudinal and transverse separation distances, respectively, and a and b are the 
coefficients of exponential decay. 


10 



The coefficient of exponential decay is a function of the particular environm rot under study 

and may also depend on free-stream Mach number as well as local flow conditions. Derivations 

of the coefficient of exponential decay are presented in Reference 26 for each unsteady flow 

condition where the data was sufficiently well defined to merit an analysis. For attached 

boundary layer flow, the coefficient of exponential decay in the longitudinal direction was 

found to be a = 0.10. For two-dimensional and axisymmetric separated flows, the values 

ranged from 0.13 at M = 2.5 to 0.33 at M =1 .60. The coefficient of exponential decay 

ao ao 

in the inner and outer regions of three-dimensional separated flows upstream of cylindrical 
protuberances were approximately the same in the longitudinal direction with a value of a ^0.7. 
in the protuberance wake, the $; tial decays were clearly defined only in the neck and far-wake 
regions with a value of a =s0.l9. Typical co- and quad-spectra are presented in Figures 8 and 
9 for attached flow and separated flow, respectively. 
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4.0 FLUCTUATING PRESSURE ENVIRONMENTS FOR SPACE SHUTTLE 


4.1 Description of Space Shuttle Vehicles 

Space Shuttle studies over the past three years have considered a wide variety of system 
confi jurations. The final concept will consist of on Orbiter vehicle with an external, 
expe ndable, liquid hydrogen/liquid oxygen system (H/O) tank. Two solid rocket boosters 
(S Sf . will be used for the shuttle booster in a strap-on configuration similar to the Titan III 
co ‘:,ept. Selection of solid rocket diameter is still under study with a 144 inch diameter 
seg Tented motor or a 142 inch diameter monolithic motor being strong contenders. Orbiter 
de ian is still under study, with the most likely candidate being a lightweight double delta 
configuration with a target gross weight of 150,000 lb. 

Paw* r for the Orbiter will be supplied by three high-pressure engines which will burn in 
pars lei with the solid rocket motors during launch. One concept for the Space Shuttle is 
t'ni aseline configuration shown in Figure 10. The detail geometry of the mated Space 
Shul'le configuration has certainly not been finalized; however, the concept has progressed 
to the point that reasonable estimates can be made in regard to the fluctuating pressure 
environments which may be anticipated during launch. The baseline configuration shown in 
rigur's 10, consisting of a lightweight Orbiter and 142 inch diameter monolithic SRB's will be 
used in this report as a baseline model for arriving at preliminary estimates of the Space 
Shuttle fluctuating pressure environments. 

4.2 Description of Flow Fielcfe 

Tf e Space Shuttle launch configuration is a complex combination of several bodies — with 
each major body having a different geometry. The flow field for the mated launch vehicle 
wi I consist of the individual flow fields of each major body of the vehicle with modified 
regions as generat- j by the interference and interactions among the flow fields of the various 
bodies. An extensive experimental/theoretical study will be required to accurately define the 
resulting flow over the mated configuration. The problem at the present time is to arrive at 
some feasible means of defining the aerodynamic flow field so that preliminary estimates of 
surface fluctuating pressure environments can be made. The approach *o be taken is to analyze, 
first, the interfe.ence-free flow field for each major body of the mated Space Shuttle vehicle, 
and to modify the individual flow fields according to predicted flow interference effects. 

Only a cursory desrriplion of the flow fields for the mated configuration is thus obtained; 
however, rhe analys. - does point out representative fluctuation pressure environments which 
may be onticipc d during Space Shuttle flight. 

The fluctr ring pressure environments are, of course, a function of launch trajectory. For the 
purpose . * this report, predictions will be made for a due east trajectory. Flight parameters 
for thw trajectory are shown in Figure 11 . It should be noted that the basic fluctuating 
er onments are functions of Mach number and scale with dynamic pressure; therefore, the 
present predictions c^n generally be corrected for other trajectories. An important exception 
to this scaling rule.* is plume induced separation, where the geometry of the separated flow, 
as well as ove'Jil level depends on dynamic pressure. Plume effects must therefore be calcu- 
lated for cacn specific trajectory. 
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4*2.1 Interference-Free Flow Environments — The solid rocket motors are basically 
cone-cylinder missile configurations withTa 17.5 degree cone half angle. Cone-cylinder 
geometries have been studied e; tenslvely in past years and their unsteady flow fields are 
fairly well defined. For these bodies, the following environments may be anticipated. 

• Attached turbulent boundary layer flow 

• Shoulder-induced separated flow (subsonic Mach numbers) 

• Terminal shock-wave oscillation (transonic Moch numbers below 1 .0) 

• Exhaust plume-induced flow separation with shock -wave oscillation 
for high altitude supersonic flight. 

For a cone-cylinder body, the types of unsteady flow conditions which will occur at a given 
Mach number depend on cone -angle, angle of attack, and Reynolds number. Typical flow- 
fields for cone-cylinders are shown for the zero angle of attack, interference -free case in 
Figure 1 2. 

For all Mach numbers, the flow over the nose cones will be attached with the location of 
boundary layer transition occurring as a function of Reynolds number. Flow characteristics in 
the vicinity of the cone-cylinder juncture will vary depending on free-stream Mach number, 
angle of attack, and shoulder ongle. For cones having half angles (shoulder angles) greater 
than approximately 10 degrees, the flow will not fully expand over the shoulder at subsonic 
speeds but will separate at the shoulder. A separation bubble is thus formed, with reattach- 
ment occurring o short distance aft of the shoulder with the axial extent of the bubble depending 
primarily on cone-angle. At high subsonic Mach numbers, the flow negotiates the shoulder 
without separating, reaches supersonic speed immediately aft of the shoulder, and produces a 
near-normal, terminal -shock wave at some point along the cylinder. The boundary layer 
becomes extremely thin at the shoulder due to the strong expansion fan and a new boundary 
layer is essentially formed at the shoulder. The boundary layer remains attached between the 
shoulder and the shock wave. At the shock location, the boundary layer may undergo separation 
depending on the strength of the shock wave. The shock strength is at o maximum corresponding 
to the initial attachment of the flow at the shoulder and decreases in strength as it moves aft 
with increasing Mach number. At sonic speed, the terminal -shock wave will move aft of the 
body and will be reduced to zero strength. 

For Mach numbers above 1 .0, the interference-free flow fields for the SRB will consist of 
attiched boundary layer flow and possible shock-oscillation environments corresponding to 
surface protuberances and other geometric irregularities. At supersonic, high altitude, 
conditions, spreading of the rocket exhaust plume will induce separated flow over the aft 
portion of the total Space Shuttle vehicle. » his will be discussed in Section 4.2.3, Plume- 
Induced Sepaiated Flow Environments. 
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H/O Tank — The H/O tank is basically an ogive-cylinder configuration. 

The configuration shown in Figure 10 has a sphere-cylinder at the nose housing a 
de-orbit engine. Design of the H/O tank is still in a fluid state/ with such details 
as de-orbit engine location and precise nose geometry not finalized at the present 
time. Latest concepts call for de-orbit engines to be located flush along the sides 
of the H/O tank, rather than in the sphere-cylinder as shown. Therefore, for the 
purposes of this report the H/O tank will be treated as an ogive-cylinder. 

Ogive-cylinders have been studied in the past, although not as extensively as the 
more basic cone-cylinder configuration. A key factor in this type of body is the 
configuration of the shoulder. 

Robertson (Reference 30) investigated the effect of introducing a radius at the shoulder 
of a cone-cylinder, with geometry ranging from a sharp corner to an ogive-cylinder. 

It was found that flow separation decreases as shoulder radius increases, with no 
separation observed for ogive-cylinders. Static pressure measurements indicate that 
the terminal shock wave on an ogive -cylinder is generally weaker and occurs at a 
higher Mach number than on an equivalent cone -cylinder. Interference -free flow 
environments for the H/O tank will therefore consist of attached boundary layer 
flow and terminal shock wave oscillation (though not as severe as for a cone-cylinder) 
at subsonic speeds, and attached boundary layer flow for Mach numbers above 1.0. 

As with the SRB, there may also be separated flow and shock oscillation environments 
corresponding to surface protuberances and other irregularities. 

Orb? ter Vehicle — The Space Shuttle Orbiter is a fairly clean vehicle from an 
aerodynamic viewpoint; however, there are several specific areas of concern which 
should be noted as regions of potentially severe fluctuating pressure environments. 

In general, a major portion of the upper surface area of the Orbiter will experience 
attached flow with relatively small fluctuating pressure levels. The lower surface 
area will experience attached flow also; however, the low^r surface will be affected 
largely by interference flow due to the close proximity of the SRB and H/O Tank. 

The following exceptions to attached flow over the upper surface c/ the Orbiter 
should be noted: 

1 ) Separated flow and shock-wave oscillation in the vicinity 
of the crew compartment. 

2) Terminal shock-wave oscillation aft of the crew compartment 
over the upper surface of the fuselage. 

3) Protuberance-induced separated flow and shock-wave oscillation 
in the vicinity of the OMS Engines. 

4) Exhaust plume-induced separated flow and shock-wave 
oscillation at high altitude supersonic flight. 
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A detailed analysis of the crew compartment and protuberance-induced fluctuating pressure 
environments will not be presented in this report since the final design of the Orbiter will 
likely be altered from the design of these components as shown in Figure 10. However, 
certair general features of these components may be expected in the final design. 

Crew-Compartment — Due to the location of the crew compartment and the require- 
ments for pilot visibility, the flow over the upper fuselage will be perturbed to a 
significant degree. The angle between the pilot's window and fuselage will form a 
compression corner and the flow in this area v/il! resemble a forward facing step or 
wedge characterized by compression-induced separated flow. At supersonic Mach 
numbers, shock-wave oscillation may be expected at the flow separation and reattach- 
ment points. Also, if the shoulder above the pilot's station is a sharp angle, over- 
expansion at transonic Mach numbers similar to that observed at cone-cylinder junctures 
may be expected. Finally, terminal shock-wave oscillation will occur over the upper 
surface of the fuselage at high subsonic Mach numbers for some distance aft of the 
crew compartment. For the baseline configuration, the shoulder above the crew 
compartment is rounded and shoulder-induced separated flow is not anticipated. 

However, terminal shock -wave oscillation will occur but with redo;:ed fluctuating 
pressure levels as compared to a sharp-shoulder configuration (see Reference 30). 

QMS Pods — The OMS units are three-dimensional protuberances that will induce 
separated flow. Relatively severe fluctuating pressures are anticipated in the vicinity 
of these units. Most of the aft fuselage and wing root section on the upper surface 
will be affected by the OMS engines. Essentially, these engines are small strap-on 
rockets and many of the arguments presented in the next section for the SRB interference- 
induced flow field are pertinent to the OMS engine areas. 

The flow field for the OMS engines will be characterized by compression-induced 
separated flow and shock -wave oscillation. 

4.2.2 Mated Vehicle Interference-Flow Environments — The launch configuration of 

Space Shuttle vehicles will be comprised of four approximately parallel bodies. Such vehicles 

have regions of flow interference between the various bodies which will result in relatively 

severe fluctuating pressure environments. For other similar configurations such as Tiran HI, 

shock-induced pressure fluctuations have been observed at the location of shock impingement 

from one body onto the surface of the ether holies. This phenomenon may be anticipated for 

Mach numbers starting at approximately 1 .2 when the bow shock for the mated vehicle divides 

into separate bow shocks for each major body of the vehicle. Shock impingement will continue 

to produce severe fluctuating pressure levels until booster separation, H/O Tank separation, 

or until the q decreases to a level sufficiently low that the disturbances are no longer of 
oo 

concern. The following shock interference regions will occur: 

1) Bow-shock impingement and shock-boundary layer interaction between 
the SRB and H/O Tank due to SRB bow shock. 
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2) Bow-shock Impingement and shock -boundary layer Interaction between 
the H/O Tank and Orblter Fuselage due to Orbiter bow shock. 

3) Shock -boundary layer interaction between the Orbiter lower wing 
surfaces and SRB due to wing leading edge shocks. 

4) Possible transonic shock-boundary layer interaction between SRB and 
H/O Tank and H/0 Tank and Orbiter due to transonic terminal 
shock waves . 

The most severe interference flow regions are anticipated to be in the vicinity of the cone- 
cylinder juncture of the SRB, in the region between the SRB and H/O Tank, and in the 
vicinity of the H/O Tank and the underside of the Orbiter nose. The most critical Mach 
number range is anticipated to be in the transonic region (0.80 < M £1.6) corresponding 
to maximum dynamic pressure, transonic buffet phenomena and near-normal shock wove 
impingement. Downstream of the shock-interaction regions, the flow may become choked 
and a shock-free interference flow environment will occur. This type environment will be 
referred to as shock -free interference flow and fluctuating-pressure levels on the order of 
those experienced in homogeneous separated flow may be expected. 

4.2.3 Plume-Induced Separated Flow Environments - - At high altitudes, a rocket 
exhaust plume expands to form an effective flared afterbody larger than the vehicle itself. 
Under supersonic conditions, flow over the vehicle may be separated due to this effective 
flare. The onset of plume-induced separation, and the extent of the separated region, is 
highly dependent on vehicle design parameters and flight conditions. Plume diameter 
increases with diminishing dynamic pressure; the onset of plume-induced seoaration therefore 
will generally occur at times after conditions of maximum dynamic pressure. As altitude 
increases after the onset of separation, the separation region grows. Pressure fluctuations 
within the separated region are similar to those for an equivalent rigid flare; howe''e r , the 
non-rigid nature of the plume boundary leads to some important differences. The most 
important difference due to motion of the plume interface is an increase in shock oscillation 
distance and attendant fluctuating pressure level. Effects of plume-induced separation have 
been investigated as part of this program. Prediction schemes for both separation length and 
fluctuating environments on axisymmetric vehicles are presented in Reference 29. 
Plume-induced separation has been computed for axisymmetric representations of the 
SRB and of the Orbiter. Tables 1 and 2 summarize vehicle and engine parameters 
pertinent to plume-induced separation on the baseline shuttle configuration. Approxi- 
mations involved are similar to those made in Reference 29 for the PRR configuration. 

Figure 13 shows the predicted separation point on the shuttle SRB and Orbiter. Separation on 
the SRB begins shortly after maximum q , and extends over most of the SRB by SRB cut-off. 

Separation on the Orbiter does not occur until well post maximum q , so that separation due 
to the Orbiter engine plume will not be an important environment, 00 Separation due to the 
Orbiter engines occurs at a much later time than that due to the SRB 's because the Orbiter's 
high H/O engines are much less underexpanded (through both higher area ratio and higher 

exhaust gamma) than the SRB solid fuel engines. 
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TABLE 1 . SRB PARAMETERS 


Vehicle Radius 

71 Inches 

Thrust 

2.5 x 10 6 Lb 

V 

' ex 

1.145 

Nozzle Area Ratio 

7.0 

Exhaust Mach Number 

2.9 

Chamber Pressure 

730 psia 


TABLE 2. ORBi i ER PARAMETERS 


Vehicle Radius 

131 Inches 

Thrust 

1.41 x 10 6 Lb 

r 

1.23 

ex 


Nozzle Area Ratio 

35 

Exhaust Mach Number 

4.05 

Chamber Pressure 

3000 psia 
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4.3 Predicted Fluctuating Pressure Levels 

The following summary presents estimates of overall fluctuating pressure levels, as a function 
of free-stream Mach number and dynamic pressure, typical of those which may be anticipated 
for the various unsteady flow environments for the Space Shuttle bodies. 

4,3.1 Attached Turbulent Boundary Layer Flow — 
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Protuberance-Flow Fields — All Mach numbers 


Interference-Flow Fields — All Mach numbers 


Figure 14 shows the basic root mean square boundary layer fluctuating pressure, both in pounds 
per square and dB, for the launch trajectory shown in Figure 11 . Maximum level for 
this environment is 1 37.5 dB at TALO = 50 sec, for which M =1.05 and altitude is 
24,400 ft. Note that this occurs before conditions of maximum Q, at 60 seconds. 


4.3.2 Shoulder-Induced Separated Flow — For the shoulder-induced separated flow 
environment, there are two sources o? fluctuating pressures which should be noted. First, 
immediately aft of the shoulder and within the separation bubble, the pressure fluctuations 
appear to be fairly homogeneous with levels similar to those observed immediately downstream 
of a rearward facing step or in the wake of bulbous payload configurations. However, further 
aft, the flow reattaches and higher level fluctuating pressures have been observed in the 
vicinity of the reattachment point. Also, even though the flow separates at the shoulder, 
the flow may reach local supersonic conditions as it expands over the separation bubble, in 
which case, reattachment aft of the shoulder will produce a recompression shock wave and 
relatively severe fluctuating pressure levels. Measurements presented in References 18 and 31 

indicate fluctuating pressure levels ranging up to yf p^ /q - 0.16 for reattachment aft 

oo 

of cone-cylinder junctures with shoulder angles ranging up to 30 degrees. Most of the results 
presented in References 18 and 31 show considerable scatter with variations in both cone 
angle and free-stream Mach number. A prediction curve which appears to be representative 
of the maximum fluctuating pressure levels for shoulder-induced separated flow on cone-cylinder 
bodies is given in Figure 15. This curve is attributed to Stevens in Reference 32, and it will be 
employed in the present predictions for shoulder-induced separated flow. Results from References 
18 and 31 are also shown in Figure 15. 


For the Mach number range from 0.60 to the attachment Mach number (M^), the homogeneous 
region of separated flow will produce fluctuating pressure levels in the range from 


P " / V 

CD 


0.02 to 0.04 (Figure 2). The fluctuating pressure levels given by Steven's prediction are 
greater than those for homogeneous separated flow and the difference can be attributed to the 
disturbances in the reattachment region (see Figure 12). For cone -cylinders with shoulder 
angles less than approximately 10 degrees, shoulder separation does not occur and the fluctuating 
pressure levels will decrease to the levels observed in attached flow as indicated by the dashed 
portion rf the prediction curve. 



Figure 15 


H/O Tank 




None anticipated for the ogi ve-cylinder configuration. 
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SRB 


• Cylinder 


Orb iter 

• Crew Compartment — Shoulder angle unknown. - anticipated for 

baseline configuration with ; d shoulder. 

Expansion-induced separated environments occur over a limited transonic Me mber range. 
Basic fluctuating levels for the Shuttle are discussed in Section 5.0. 

4.3.3 Transonic Shock-Wave Oscillation — Typical cone-cylinder data showing the 
attachment Mach number and the variation of shock location with free-stream Mach number is 
presented in Figure 16. Shown in Figure 17 is an interpolation of these data for two additional 
core angles. Shock-induced fluctuating pressure levels, normalized by free-stream dynamic 
pressure, are presented in Figure 18 for a range of cone angles. Curves for 17.5° and 22.5° 
are interpolated, 17.5° corresponding to the SRB. At this time, corresponding predictions 
have not been made for the H/O tank. Fluctuating levels and X^/D fo- - an or. ve-cy Under 

are generally less than for a cone-cylinder of similar nose angle; however, geometry of the 
H/O tank has not yet been finalized to a point where precise predictions are reasonable. 

The nose half angle of the ogive in Figure 10 is approximately 38° . As a preliminary estimate, 
fluctuating levels will be taken as corresponding to a 30° cane. This is a very conservative 
estimate, and will verpredict the level by at least several dB. X /D is not estimated for 
the H/O Tank. SW 

Transonic shock-wave oscillation will also occur aft of the Orbiter crew compartment canopy 
with characteristics similar to H/O tank and SRB. Because of unknowns associated with canopy 
geometry, precise estimates of fluctuating pressure levels and shock-wave locations as a function 
of Mach numbers are not made herein. However, the data for the H/O tank and SRB are 
indicative of the ranges expected for the Orbiter. 

The prediction for shock oscillation is. 

/ a a, ' f(M oo' 0 N ) ' F^re 17 
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H/0 Tank 


• Cylinder 
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• Cylinder 


Orb iter 
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in Figure 16. 
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approximately 3 6^ width 

correspondin'- to X /D 
sw 

in Figure 17. 


Location and Mach number 0.86 M <1.0 
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approximations of the range 
of parameters for Orbiter 
terminal shock-wave 
oscillation. 


Transonic shock oscillation occurs over a limited Mach number range. Fluctuating levels and 
shock location are easily obtained from Figures 16 through 18; levels and location for the 
shuttle are computed and discussed in Section 5.0. 

4.3.4 Compression-Induced Separated Flow (Two-Dimensional) 
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Obiter — Crew Compartment Canopy ) 

Protuberance -Flow Fields j All Mach Numbers 

Interference-Flow Fields ) 

Fluctuating levels for the launch trajectory are shown in Figure 19. Note that because 
^ P 2 /q is an increasing function of Mach number, the critical (maximum fluctuating 

level) compression- induced separation environment occurs at TALO = 72 seconds, a later time 

than maximum q . 

oo 

4,3,5 Supersonic Shock -Boundary Layer Interactions 
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Orb iter 
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Fluctuating levels for the launch trajectory are shown in Figure 20. Because the prediction 

used here is not a function of Mach number (other than applying only for M > 1), the critical 

shock oscillation environment occurs at maximum q . 

oo 

4.3.6 Plume-Induced Separation — Plume-induced separated flow is no! a basic 
environment as those Discussed in Sections 4.3.1 through 4.3.5, but is rather a modification 
of compression corner basic environments. The modification is dependent on the particular 
vehicle and trajectory. Not onlydoes the mean separated flow depend on vehicle parameters, 
but the response of the plume enters into the fluctuating environments. In Reference 29, it 
was shown that homogeneous separated flow fluctuations in plume-induced separation are 
essentially the same as for a rigid fiare, while shock oscillation depended on plume parameters. 
The homogeneous environment is thus given by Figure 19, with the (imitation that it applies 
only to flight conditions where plume-induced separated flow exists, as defined in Figure 13. 
Note that separation due to the orbiter engines occurs only at very high altitudes; the dynamic 
pressure at these altitudes is sufficiently smr hat orbiter plume-inducea separation may be 
neglected relative v o other environments. Toe SRB plumes induce separation at considerably 
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lower altitudes, and must be considered. In Reference 29, a prediction scheme was developed 
for the fluctuating pressure associated with the separation shock. Predictions made for the PRR 
configuration indicate that shock excursion and fluctuating pressure level for plume-induced 
separation is generally less than 25% greater than for a rigid flare me baseline shuttle is 
similar; shock oscillation fluctuating pressure levels for plume-induced separation on the 
baseline shuttle SRB wil! therefore be no more than 2 dB greater than for the basic shock 
environment given in Figure 20. This environment occurs only at TALO greater than 62 seconds, 
the onset of plume-induced separation. 

4.4 Predicted Power Spectra 

The following empirical prediction equation, presented in normalized form, may be employed 
to compute the power spectra for the various fluctuating pressure environments. 

4.4.1 Attached Turbulent Boundary Layer Flow — The following equations apply to 
the homogeneous regions of attached flow. 
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x = Downstream distance from the leading edge 


R = Reynolds number = U x/v 
t co 


v = Kinematic viscosity 

4. 4.2 Compression/Expansion-Induced Separatee Flow and Shock-Free Interference 
Flow — The following equations apply to the homogeneous regions of separated flow. 
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, for expansion induced 
separated flows. 


- 

the results as determined in Figure 2 for compression 
induced separated flows. 


and the subscripts l and ao refer to local and free-stream conditions, respectively. 

It is anticipated that the above equation can be used with good uccuracy to predict the power 
spectra for fluctuating pressures within the homogeneous region of expansion-induced separated 
flows and shock-free interference flows although it was derived based on data taken in 
compression corners. 
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4.4.3 Shock-Osci I lotion Environments — To predict the power soectra for non- 
homogeneous flows generated b)T shock -wave oscillation, prior knowledge of the overall 
fluctuating pressure levels is required. Under the assumption of statistical independence 
between the various contributing sources, the power spectra for nonhomogeneous environments 
may be written as the summation of power spectra of the contributing sources. 
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The power spectra for nonhomogeneous attached flow as caused by shock -wave oscillation in 
the vicinity of an aftached flow region has the following form normalized by local conditions 
upstream of the shoe!: wave: 



The following expression has been derived for nonhomogeneous separated flow: 
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also the subscripts and superscripts denote the following: 


Subscripts: SW 

S 
A 

Superscripts: 1 

H 


shock wave 
separated flow 
attached flow 

absence of viscosity (inviscid) 
homogeneous flow 


Generally, shock-wave oscillation occurs in the presence of separated flow for critical 
environments, and examples are shock-induced interference flows, terminal shock-wave 
oscillation, the separation and reattachment points for compression corner flow at supersonic 
Mach numbers, and the reattachment point for expansion shoulder flows for local supersonic 
conditions. For terminal shock-wave oscillations at Mach numbers near 1 .0, the flow does 
not separate and the nonhomogeneous attached flow equation is applicable. 
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5.0 ANALYSIS OF CRITICAL SPACE SHUTTLE FLIGHT CONDITIONS 


The analysis presented In Section 4.0 revealed that many severe and complex fluctuating 

pressure environments may experienced by the Space Shuttle vehicles. The most severe 

environments will occur at transonic and low supersonic Mach numbers, 0.60 < M < 2.0. 

ao 

In terms of fluctuating pressure coefficient, transonic environments (0.60 < M < 1 .0) are 

ao 

the most severe, while supersonic environments are critical because of high dynamic pressure. 
Although critical (maximum) fluctuating levels for each supersonic environment do not 
necessarily occur at conditions of maximum dynamic pressure (MAX Q) the levels at MAX Q 
are sufficiently close to maximum that for the purposes of this report critical supersonic 
environments have been computed at MAX Q. Transonic environments have been computed 
over a range of Mach numbers. 

The various prediction schemes developed under this contract have been incorporated into 
several computer programs. Descriprions of the various programs are contained in Appendix 
D. All predictions contained in this report have been made utilizing these programs. 
Calculations are presented for overall levels and typical spectra. Although a computer program 
has been prepared for the calculation of correlation functions (see Appendix D for description), 
specific sample calculations are not presented. The scaled correlations shown in Figures 8 and 
9 adequately illustrate these for the purposes of the present report. 

Computations were performed for the following critical flight conditions: 

1 . Maximum Dynamic Pressure (MAX Q), M = 1.31 

ao 

• Attached Turbulent Boundary Layers 

— Typical H/O Tank, SRB and Obiter Stations based on 
representative values of local boundary layer thickness 
and velocity. 

• Homogeneous Compression-Induced Separated Flows 

— Typical of the shock-free region of separated flow over the 
Crew Compartment and shock-free Interference flow fields. 

• Shock -Bound ary Layer Interactions 

— Typical of Boundary Layer Separation and reattachment 
areas and shock-induced interference flow fields. 

• Protuberance Induced Flows 

Typical of various H/O Tank, SRB and Orbiter 

protuberances 
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2 . 


Critical Expans ion -Induced Separated Flows 


• SRB Worst -Case Shoulder Separation ~ M = 0.75 

CO 

3. Critical Transonic Shock Oscillation Environments 

• SRB Cylinder - M -0.78, 0.82, 0.86, 0.90, 0.94 

oo 

• H/O Tank Cylinder - M = 0.90 

00 

• Orbiter Upper Fuselage (in the range of SRB and H/O Tank data) 

4. Critical Interference-Flow Environments 

• Typical SRB, H/O Tank and Orbiter Stations 

— Shock-Induced Interference Flow - M = 1 .31, 1 .6, 2.0, 3.0 

ao 

— Shock-Free Interference Flew - M = 0.60, 1 .00, 1 .31 

ao 

Primary emphasis has been given to the condition of maximum dynamic pressure since this flight 

condition will generally produce the most severe fluctuating pressure environments. However, 

it is equally important to consider certain critical flow conditions which are unique to Mach 

numbers below the condition of MAX Q, such as transonic shouider-induced flow separation 

and transonic shock-wave oscillation. Finally, critical interference-flow environments will be 

generated due to bow-shock impingement which will occur at supersonic Mach numbers. For 

the latter case, a range of supersonic Mach numbers are considered, starting at M =1 .31 

(MAX Q)and ranging up to M - 3.0. 00 

oo 

The computations represent typical overall fluctuating pressure levels and one-third octave band 
spectra which may be anticipated within the perturbed flow region under analysis. It should be 
noted that the predictions are based on estimated average values of local flow parameters rather 
than actual measured values. The greatest uncertainty in the predictions is in the frequency 
scale rather than the amplitude scale since errors in the specification of local flow conditions 
(typical length scale and velocity scale) will cause a translation in frequency rather than 
amplitude. If actual local values of boundary layer thickness, flow separation length and 
velocity become available, the frequency scale for the predictions may be adjusted accordingly. 
The following relationship will apply for correcting the frequency scale of the predictions. 

f = T~ TT # f P ' Hz 

p r 


where the subscript p denotes the values used in the prediction. 
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Variations in overall fluctuating pressure levels and one-third oc*ave band spectrum levels may 
be anticipated for full-scale dynamic pressures different from those employed in the predictions. 
An adjustment to the predicted values for different values of q can be made as follows: 


OAF PL 


(OAFPL) + 20 log 
P 10 



dB 


1/3 OBFPL = 


(1/3 OBFPL) + 20 log 
p to 



, dB 


In the following sections, the computed overall fluctuating pressure levels and one-third 
octave band spectra are discussed for the various fluctuating pressure environments which are 
expected to be encountered at critical Space Shuttle flight conditions. 

5.1 Maximum Dynamic Pressure, MAX Q 

The following trajectory parameters were employed in the predictions of MAX Q fluctuating 
pressure levels and spectra. 

• Mach number. M = 1 .31 

oo 

• Dynamic Pressure, q =649.4 lb/ft 2 

oo 

• Altitude, h = 34,400 ft 

• Free-Stream Velocity, U = 1,296 ft/sec 

00 

Representative spectra for the unsteady flow fields which will be encountered by the mated 
Space Shuttle vehicles at MAX Q are presented in Figure 21 . The following environments 
were analyzed using a nominal value of local boundary layer thickness of 1 .0 foot and a 
nominal value of local velocity equal to free -stream velocity. 

• Attached Turbulent Boundary Layer Flow 

• Homogeneous Compression-Induced Separated Flow 

• Shock-Boundary Layer Interactions 

• Three-Dimensional Protuberance Inc :ed Flow 
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Attached turbulent boundary layer flow may be anticipated for much of the vehicle surface at 
MAX Q with overall fluctuating pressure levels near 137 dB. With the exception of inter- 
ference flow regions between the various bodies, local protuberance flows, and exhaust plume- 
induced separated flow, all surfaces will experience fluctuating pressures characteristic of the 
attached flow environment. 

Homogeneous separated flow regions (separated flow regions not affected by shock-wave 
oscillation) will experience fluctuating pressure levels near 148 dB. Typical areas are the 
interference flow region between the H/O Tank and SRB cylinders downstr am of the shoulders, 
and between the H/O Tank and Orbiter fuselage well downstream of the Orbiter nose. Also, 
local surface irregularities such as the crew compartment canopy and two-dimensional type pro- 
tuberances will have regions of homogeneous separated flow. Shock-boundary layer interaction 
regions will experience fluctuating pressure levels near 161 dB at MAX Q and typical areas 
are the interference regions in the vicinity of the nose cones and cone-cylinder shoulders 
on the SRB and the interference region between the Orbiter nose and tj/O Tank. 

Also, boundary layer separation and reattachment points in a local supersonic environment, 
such as the crew compartment canopy and various protuberances, will experience shock-induced 
fluctuating pressures near 160 dB. Three-dimensional protuberances are defined as protuberances 
which will cause the perturbed flow to move around the protuberance transverse to the local 
unperturbed flow direction. Extremely large fluctuating pressures may be anticipated for these 
environments due to the large vortex action within the three-dimensional separated flow field 
(see Reference 22). Levels at MAX Q may be anticipated near 168 dB. 

The frequency charac .ristics of the various fluctuating pressure environments are shown in 
figure 21 and may be summarized as follows: 

• Attached Turbulent Boundary Layers — HIGH FREQUENCIES 

• Homogeneous Separated Flows (2-D and 3-D) — INTERMEDIATE FREQUENCIES 

• Shock-Boundary Layer Interactions — LOW FREQUENCIES 

Naturally, consideration must be given to local length and velocity scales when determining 
spectra for specific flow conditions. A composite of predicted spectra for various specific areas 
of the Space Shuttle is presented in Figure 19 of Reference 33 where consideration was given 
to anticipated local flow conditions and typical length and velocity scales; the present 
predictions are similar. 
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5.2 Critical Expansion-Induced Separated Flows 

Expansion-induced separated flows will occur at the cone-cylinder junctures of the SRB, 
possibly above the crew compartment of the Orbiter (depending on shoulder geometry) and 
In the wake of certain protuberances and other surface irregularities. With the exception 
of bluff geometries, expansion-induced separated flow will be confined, primarily, to 
moderate and low subsonic Mach numbers. Typical examples are the cone-cylinder junctures 
of the SRB. For these cases, maximum disturbances will occur immediately prior to flow 
attachment. A typical spectrum is shown in Figure 22. Predicted levels are near 157 dB for 
the SRB. Disturbances in the vicinity of flow reattachment aft of the shoulder peak in the 
frequency range from 8 to 1 6 Hz . 

5.3 Critical Transonic Shock-Oscillation Environments 

Following flow attachment at the shoulders of the SRB and aft of the crew compartment on the 
Orbiter, local supersonic flow will result in a near normal, transonic, terminal shock wave. 
Typical spectra for the SRB are shown in Figure 23. Flow attachment is predicted at M -0.78 

QD 

for the SRB shoulder and corresponds to the condition of maximum transonic shock-induced 

pertubations with an overall level of 164 dB. As Mach number is increased, the shock wave 

will move aft on the SRB cylinders and decrease in strength with a predicted level of 143 dB 

at M - 0.94. The frequency of peak fluctuating pressure level also decreases as the shock 
00 

moves aft due to: 1) the Increase in inflowing boundary layer thickness upstream of the shock 

wave, and 2) the decrease in local velocity which results from the isentropic recompression 

between the shoulder and the shock wave. Peak frequencies of shock-induced disturbances 

are estimated to be 40 Hz at M - 0.78 f decreasing to 4 Hz at M = 0.94, As Mach 

00 00 

numbers approach 1 .0, the shock is further weakened with surface disturbances approaching 

those of attached flow. Based on a 30° cone, as discussed earlier, critical terminal shock 

oscillation on the ogive-cylinder H/O tank is estimated to be 156.6 dB at M - 0.90. This 

oo 

level should be regarded as an upper bound, with the actual level several dB lower. 

Without detailed knowledge of the Orbiter crew compartment geometry, it is difficult to 
determine the exact nature of terminal shock -oscillation on the Orbiter fuselage. However, 
it will occur and may be expected to exhibit characteristics similar to the H/O Tank and 
SRB environments. 


5.4 Critical Interference-Flow Environments 

It was previously noted in Section 4.0 that several regions of the mated Space Shuttle vehicles 
will experience interference flow environments caused by the close proximity of the various 
bodies composing the launch configuration. Interference-flow fields on the Space Shuttle will 
be caused, in general, by the impingement of shock waves from one body onto the surface of a 
neighboring body. Naturally, shock-wave impingement is possible only in a local supersonic 
flow environment. At local subsonic Mach numbers, say downstream of shock -boundary layer 
interaction regions for supersonic free-stream conditions, the perturbed flow field will not 
contain shock waves even though the flow field is caused by interference betv;een flow fields. 
Thus, it is convenient to define two basic types of interference flow fields: 
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• Shock-Induced Interference Flows 

• Shock-Free Interference Flows 

With the exception of possible transonic terminal shock-wave impingement* , shock-induced 
interference flows will occur only at supersonic free-stream conditions whereas shock-free 
interference flows may occur at all Mach numbers. Typical spectra for shock-induced 
interference flows typical of the interference regions between the H/O Tank and SRB (due to 
SRB bow shock impingement) between H/O Tank and Orbiter fuselage (due to Orbiter bow 
shock impingement) and between the SRB and the lower surface of the Orbiter wings (due to 
the impingement of the Orbiter wing leading edge shocks) are shown for a range of supersonic 
Mach numbers in Figure 24. Also shown in Figure 24 are representative spectra for shock-free 
regions of the interference-flow field areas which are predicted to be affected by this 
environment and are the neighboring surfaces cf the SRB and H/O Tank cylinders (downstream 
of the shoulders), the neighboring surfaces between the H/O Tank and Orbiter well downstream 
of the Orbiter nose, and the lower surface of the Orbiter wings and aft cylinders of the SRB, 

5.5 Comparison of Predictions with Experimental Measurements 

In Reference 33, a detailed comparison was made between predictions for the MSC -049 Shuttle 
configuration with measurements reported in Reference 34. The same general discussion applies 
to the present predictions, and will not be repeated herein. Shown in Figure 25 is a comparison 
of present predicted levels to measured levels at critical transonic conditions. 


Recently published data (Reference 34, discussed in Reference 33) indicates that 
shock-induced interference flow due to a blocked terminal -shock wove between the 
H/O Tank and O 'biter nose does occur at near '-sonic Mach numbers. 
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6.0 CONCLUDING REMARKS 


Preliminary estimates of Space Shuffle flucfuafing pressure environments have been made 
based on prediction fechniques developed by Wyle Laboratories. Particular emphasis has 
been given to the transonic speed regime during launch of a paraMel-burn Space Shuttle 
configuration. The baseline configuration, together with a typical flight trajectory, have 
been used as models for the predictions. Analyses have consisted of the following: 

• Description of Flov Fields and Attendant Fluctuating Pressure Environments 
for the Space Shuttle Launch Configuration . In most cases, these analyses 
are sufficiently general to be applicable ro any parallel-burn launch 
configuration consisting of a reusable Orbiter vehicle, external H/O Tank 
and two strap-on solid rocket motors. 

• Specification of Interference-Free and Interference-Induced Flow Fields 
and Empirically Determined Attendant RMS Fluctuating Pressure Levels and 
Power Spectral Densities. Particular emphasis has been given to the transonic 
Mach number regime during launch. 

• Prediction and Evaluation of Plume-Induced Separated Flow Environments. 

• Analyses cf Overall Fluctuating Pressure Levels (Expressed in Decibels) and 
One-Third Octave Band Spectra of Various Fluctaufing Pressure Environments 
which are Expected to be Encountered at Critical Space Shuttle Flight 
Conditions. 

• Comparison of Analytical Predictions with Wind Tunnel Measurements for 
Various Regions of a Typical Space Shuttle Vehicle. 

The results of this study have resulted in the following conclusions* 

1) The following critical flow environments together with typical overall fluctuating 

pressure levels are predicted: 

• Cone-Cylinder Shoulder-Induced Separated Flow on the SRB with OAFPL 
Peaks Near 157 dB for Zero Angle of Attack 

• Terminal Shock-Wave Oscillation on the SRB, H/O Tank and Orbiter 
Fuselage Upper Surface Aft of Crew Compartment with OAFPL Peaks 
Near 164 dB 

• Three-Dimensional Protuberance Flows Representative of Various Space 
Shuttle Protuberance Environments with OAFPL Peaks Near 168 dB 

• Compression-Corner Induced Separated Flow in the Vicinity of the Orbiter 
Crew Compartment with OAFPL Ranging up to 148 dB and Peaks Near Flow 
Separation and Reattachment Ranging up to 157 dB 
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• Interference-Flow Environments Between the SRB and H/O Tank, Between 
the H/O Tnnk and Orbiter Ne«e, and Between the Orbiter Wind and 

SRB with Peak oAFPL Ranging up to 161 dB 

• General Attached Flow Environments with OAFPL Runqing up to 137 dB 

2) One-third octave band spectra for the above environments were generally one of 
three types: 

• Attached Turbulent Boundary Layer Spectra — Typically High Frequencies 

• Homogeneov: Separated Flow and Shock-Free Interference Flow Spectra — 
Typically Intermediate Frequencies 

• Shock-Oscillation and Shock-Induced Interference Flow Spectra — Typically 
Low Frequencies 

3) Flow separation due io the SRB engine exhaust plumes will occur shortly after MAX Q. 
Fluctuating environments will be comparable to those associated with compression 
corners. Separation due to the Orbiter engine exhaust plume occurs at high enough 
altitude such that it does not represent an important fluctuating pressure environment. 

4) Comparisons between predicted and experimentally measured OAFPL showed good 
agreement except over aft portions of the Orbiter fuselage. In gene'al, the agreement 
between OAFPL values for experiment and predictions was within 3 dB. 
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FIGURES 





Figure 1 . 


Comparison of Pressure Fluctuation Measurements by Various Investigators 

for Attached Turbulent Boundary Layers 
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Coe and Rechtien, 45° Conical Frustrum 



40 


Figure 2. Variation of Fluctuating Pressure Level for Homogeneous Separated FI 
wifh Mach Number 



Attached Flow 



Figure 3. Variation of Fluctuating Pressure Level for Shock-Wave Oscillation with 
Mach Number 








Figure 4. Power Spectra for Turbulent Boundary Layer Fluctuating Pressures 
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Figure 7. Longitudinal Distribution of Pressure Fluctuations and Typical Power 

Spectra in Vicinity of Supersonic Flow Separation Ahead cf a 45° Wedge 
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Figure 10. Schematic of Baseline Shuttle Launch Configuration 
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Shoulder-Induced Separated Flow 
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Figure 12, Schematic of Interference-Free Flows Over 
a Cone -Cylinder Body 
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Figure 14 . Basic Attached Flow Fluctuating Pressure Level, Shuttle Launch Trajectory 
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One-Third Octave Band Fluctuating Pressure Level, dB 





One-Third Octave Band Fluctuating Pressure Level, dB 



One-Third Octave Band Center Frequency, Hz 


Figure 22. One-Third Octave Band Spectrum for Shoulder-Induced 
Separated Flow on the SRB at Low Transonic Mach 
Number 
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One-Thire Octove Band Fluctuating Pressure Level, dB 



One-Third Octave Band Center Frequency, Hz 

Figure 23. One-Third Octave Band Spectra for Terminal Shock-Wave 
Oscillation on the SRB at Transonic Mach Numbers 
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One-Third Octave Band Fluctuating Pressure Lev'el, dB 
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Figure 25. Comparison of Meosured and Predicted Overall Fluctuating Pressure 
Levels for a Typical Porollel-Bum Space Shuttle Launc Configura- 
tion at Transonic Mach Numbei 
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APPENDIX A 


SUMMARY OF PREDICTION METHODS FOR 
IN-FLIGHT FLUCTUATING PRESSURE ENVIRONMENTS 

.0 INTRODUCTION 

When a vehicle moves through air there are two basic means by which it can produce 
noise: (1) by its propulsion mechanism (motor-jet, rocket, etc.,) and (2) by its 
interaction with its surroundings. At low speeds, for example, during and immediately 
after lift-off, the first of these is by far the dominant one while near or above the 
speed of sound, mechanism (2) becomes most important. 

During any flight cycle for an aerospace vehicle, there are three important phases of 
the flight which should be investigated ir order to assess the structural loading due to 
fluctuating pressures. These are listed in the chronological order in which they occur. 

e Lift-off phase during which acoustic excitation results from the rocket 

exhaust noise. 

• Launch flight to orbit phase, during which rocket exhaust noise diminishes 
and aerodynamic fluctuating pressures (pseudo-sound) starts to dominate. 
From an aerodynamic noise viewpoint, this phase becomes most critical 

at transonic Mach numbers (0.60 £ M < 1 .6) 

• Re-entry phase during which only aerodynamic fluctuating pressures are 
present. 

This Appendix is devoted to the specificarion of surface fluctuating pressures resulting 
from unsteady aerodynamic phenomena during the launch phase of flight. Aerodynamic 
fluctuating pressures (pseudo-sound) are zero at launch and increase to peak values as 
the vehicle passes through the transonic Mach number range. Previous wind tunnel 
and flight data show that fluctuating pressures are proportional to free-stream dynamic 
pressure q (= y P M* /2 , where y is the ratio of specific heats, P is the free- 

QO 00 00 00 

stream static pressure, ar.d M^ is the free-stream Mach number) for a given unsteady 
flow h .ienomenon. However, peak fluctuating pressures do not necessarily occur at 
maximum q for certain regions of a vehicle c je to the non-homogeneous nature of 

GO 

the flow field. For example, regions of thu vehicle exposed to separated flow and the 
impingement of oscillating shock waves will experience fluctuating pressures at least an 
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order of magnitude greater than regions exposed to attached flow. Thus, if separated 

flow and oscillating shock waves are present, say at Mach numbers other than the range 

of maximum q , then peak fluctuating pressures will also be encountered at conditions 
00 

other than at maximum q^ . Thus, it is easily seen that vehicle configuration is very 
important in the specification of fluctuating pressure levels since the source phenomena 
are highly configuration dependent in addition to varying with Mach number and angle 
of attack. 

In light of the foregoing discussion, one general statement can be made in regard to 
aerodynamic fluctuating pressures. Regions exposed to the same unsteady phenomenon 
will experience fluctuating pressure levels which are proportional to free-stream 
dynamic pressure. Thus, it can be readily seen that a fundamental parameter in the 
specification of the surface excitation is free-stream dynamic pressure and its variation 
with Mach number. For a given configuration, Mach number and angle of attack define 
the phenomena, and dynamic pressure defines the fluctuating pressure levels associated 
with the phenomena. 

Unsteady aerodynamic flow and the attendant fluctuating pressures experienced by 
aerospace vehicles naturally depend on the flight environments and the geometry of 
the vehicle. There are an infinite number of possible configurations and any dis- 
cussion of their fluctuating pressure environment must be general. Practically all 
experimental data for unsteady aerodynamic flow have been : red for bodies of 

revolution which are typical of missile configurations. As a resuit of these studies, 
it is well known that certain basic unsteady flow conditions will occur regardless of 
the detailed geometry of the vehicle. The occurrence of these basic fluctuating pres- 
sure phenomena and their statistical properties can be predicted quite accurately. It is 
convenient to discuss these basic flow conditions for bodies of revolution; however, 
this is certainly no restriction on either the feasibility or the practicality of predicting 
their occurrence on more complicated configurations. Thus, in the following paragraphs, 
general features of typical bodies of revolution are defined and the unsteady flow fields 
which they encounter are discussed. Furthermore, uerospace vehicles may have a number 
of protuberances projecting from their surfo-- in which case the flow field is complicated 


A-2 



by the super-posifi^n of the protuberance flow field onto the flow field of the basic 
structure. Most protuberances are three-dimensional projections and general charac- 
teristics of these flows should be considered as separate and unique problems. 

2.0 BASIC FLUCTUATING PRESSURE PHENOMENA 

Examples of several bodies of revolution are shown in Figure 1 . For the purpose of 
the present discussion, three basic configurations will be considered as specified 
below: 

• Cone-cylinder shroud 

• Cone-cylinder-flare shroud 

• Cone-cylinder-boattail shroud 

Virtually all axisymmetric vehicles fall into one of these categories although numerous 
modifications to the basic geometry have been employed in the past. 

Several fluctuating pressure environments having different statistical properties may 
exist over a vehicle at any given instant in the flight trajectory, it is convenient to 
consider three separate Mach number ranges — subsonic, transonic, and supersonic — 
for each of the three basic shroud configurations. Further, the flow fields will depend 
on the angle of attack of the vehicle which causes nonsymmetrical loading (both 
statically and dynamically); however, for the purpose of this discussion, nonsymmetri- 
cal loading will not be discussed. 

Schematics of subsonic, transonic and supersonic flow fields for the basic configurations 
are shown in Figure 1 . At subsonic speeds, all three configurations experience regions 
of attached flow and separated flow. The cone-cylinder portion of each configuration 
induces separated flow immediately aft of the cone cylinder juncture for cones having 
half-angles greater than approximately 15 degrees. Re-attachment occurs within 
approximately one diameter aft of the shoulder (depending on cone angle) for the cone- 
cvlinder and boattail configurations, whereas for the flare body, separation may con- 
t.nue over the flare. Both the flare and boattail induce separation for typical con-' j - 
urations. At high transonic speeds, the flow negotiates the shoulder of a cone-cylinder 
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body without separating, reaches supersonic speed Immediately aft of the shoulder and 
produces a near-normal, terminal, shock wave a short distance aft of the shoulder. 

The boundary layer immediately aft of the shock may or may not separate depending 
on the strength of the shock wave. At transonic speeds, the boattail and flare region 
produce separated flow which may be accompanied by weak shock waves in the vicinity 
of the separation and reattachment points. At supersonic speeds, the cone-cylinder 
configurations produce regions of attac L ed flow. For the flare configuration, the 
separated flow is bounded by shock waves at the separation and reattachment points, 
whereas for the boattail configuration, separation occurs at the shoulder of the boat- 
tail (expansion region) and is bounded at the reattochment point by a shock wave. 

It is evident that even simple vehicle shapes, such as cone-cylinders, produce complex 
and highly non homogeneous flow fields at certain Mach numbers — particularly at 
subsonic and transonic speeds. The unsteady flow phenomena are of particular 
importance at transonic speeds, since in this range, fluctuating pressures r?ach maxi- 
mum values due to their proportionality to dynamic pressure. In order to assess the 
fluctuating pressure environment of a vehicle cf any arbitrary geometry, it is conven- 
ient to discuss the statistical properties of the fluctua.ing pressures for each of the basic 
types of unsteady flow condition. From Figure 1 it will be noted that the following 
flow conditions may occur for various regions of a vehicle. 

• Attached flow 

• Separated flow 

• Shock-boundary iayer interaction 

Each of the above flow conditions exhibits different statistical characteristics. 

Attached flow pressure fluctuations result from the disturbances within turbulent 
boundary layers. Separated flow pressure fluctuations result from disturbances within 
the separated shear layer and instabilities associated with the separation and reattach- 
ment points. Pressure fluctuations for shock-boundary layer interaction result from the 
movement of the shock wave and the static pressure discontinuity associated with the 
shock wave. The statistical characteristics of each fluctuating pressure environment 
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that are important in the analysis of structural response may be classified under three 
parameters: 

• The overall level 

• The power spectrum 

• The cross-power spectrum (or narrow band cross correlation) 

Each unsteady flow condition with general statistical characteristics will be discussed 
separately in the following subsections. 

2.1 Attached Turbulent Boundary Layers 

The surface fluctuating pressures beneath attached turbulent boundary layers have been 
the subject of both theoretical and experimental study for a number of years. The 
turbulent boundary layer extends over a considerable portion of the surface of vehicles 
in flight and, thus, it is considered to be one of the principle sources of aero-acoustic 
excitation to the vehicle structure. Several years ago, workers such as Kraichman, 
Lilley, and Hodgson deve loped theoretical formulations for the fluctuating pressures 
under urbulent boundary layers and, more recently, several carefully planned experi- 
ments have provided additional information on the statistical characteristics of the 
pertubations . Lawson, Reference 1 , presents a good summary of the results of studies 
on this subject, with the exception of some recent measurements by NASA-Ames. In 
Lowson's report, the basic mechanism underlying the production of the surface pressure 
fluctuations beneath turbulent boundary layers is discussed, together with a presentation 
of empirical and semi-empirical prediction techniques. This section of the present dis- 
cussion is a brief overview of Lowson's prediction formulae with the exception of the 
power spectra, which has been modified to be more consistent with the power spectra 
at low 3trouhal numbers. The following discussion presents a review of the experimen- 
tal results and prediction formulae in terms of the most important statistical parameters. 
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Overall Level 


The correct method of presenting overall fluctuating pressure levels for surfaces 
beneath the convected turbulence in boundary layers is in terms of the root-mean- 
square fluctuating pressure level. vr. Free-stream dynamic pressure, q^ , local 


ao 

dynamic pressure, q , and wall shear stress, r , have been used to normalize 

I — - l w 

V p so that meaningful data collapse can be realized throughout the Mach number 

range. The most generally accepted normalizing parameter is q^ and thus, will be 

used in the current expressions. 


The effects of free-stream Mach number, M , on the normalized RMS intensities of the 

CD 

fluctuating pressures in attached flows are shown in Figure 2. There is significant 
scatter in the data which may be attributed to several factors: 1 ) background noise and 
free-stream ft lence in the testing medium, 2) instrumentation quality and the preci- 
sion of the experimental technique, 3) data acquisition and reduction techniques, etc. 
For the range of Mach numbers covered in the data of Figure 2, the normalized RMS 
value of the fluctuating pressure varies from V P 2 / q^^ ^*006 subsonic Mach 
number* to 0.002 at supersonic Mach numbers. Lowson, Reference 1, proposed !*he 
following semi-empirical prediction formula which appears to agree with the general 
trend in the data: 


V7^/q = 0.006/(1 + 0.14 M 2 ) (1) 

/ 00 00 

If is important to note that this formula has some theoretical basis and is not strictly 
an empirical approximation of measured results (see Reference 1 ). The use of this 
formula at high supersonic and hypersonic Mach numbers should be done so with the 
understanding that it has not been verified in this Mach number range and may lead to 
significant error. However, in the Mach number range up to, say M = 3.0, it 

OD 

is in good agreement with experimental results. 



It should be noted that the results presented in Figure 2 , particularly the wind tunnel 
results, were obtained for both homogeneous and stationary flows at free-stream con- 
ditions and in the absence of external pressure gradients. Consideration should be 
given to local conditions which deviate from free-stream conditions. 


Power Spectra 

Power spectra represent the distributions of the mean square fluctuating pressure with 
frequency. Power spectra for attached turbulent boundary layers are found to scale 
on a Strouhal number basis; that is, the frequency is normalized by multiplying by a 
typical length and dividing by a typical velocity. The advantages of using normalized 
spectra are obvious since it enables similar, homogeneous, flows to be represented by 
a single spectrum regardless of the scale of the flow field or the free-stream velocity. 


Numerous studies have been conducted to determine the proper parameters to be used 
to nondimensionalize the spectra for various aero-acoustic environments. Unfortunately, 
the choice of parameters which best collapses the data appear to be dependent on the 
nature of the fluctuating pressure environment. In general, free-stream velocity is 
used as the normalizing velocity parameter, although a typical eddy convection 
velocity (itself a function of frequency) has been used occasionally. The local con- 
vection velocity appears to correspond more closely with the physical situation for 
fluctuating pressures due to turbulent eddies. Selection of a typical length is more 

difficult. Boundary layer thickness (6, ), displacement thickness (S*), wall shear 

b 

stress (t ) and momentum thickness (0) have all been used by various investigators. 

The most generally used typical lengths ere 6, and 5*. 

b 


Lawson, Reference 1 , proposed an empirical formula for the power spectrum for 
attached turbulent boundary layers based, primarily, on the experimental results of 
Speaker and Ailman. In comparing this formula with other data, and in particular, 
with recent measurements at supersonic speeds by NASA-Ames, the Lowson prediction 
appears to underestimate the spectral levels at low Strouhal numbers and also gives 
too large a roll-off at high Strouhal numbers. Therefore, a new formula is presented 
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which appears to be more representative of e <perimental findings throughout the Mach 
number range. In this formula, it will be noted that 6* and U (the free-strearr. 

CO 

velocity) have been used as normalizing parameters. The power spectral density, 

0(u) is given by the relation: 


0(u) U 

CD 




u 6* 
0 

IT - 


I 

I 


0.9 

I + ( u/ u 0 ) 


| 2.0 
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(2! 


u 

where u =0.5 — ^ 


(0.006) 2 


(1 + 0.14M 2 f 


5* = 6^/8 for M < 1 .0 


(1.3 + 0.43M 2 ) 5 

6* = — j — for M> 1.0 

10.4+ 0.5M^ [l + 2 . 1 0“ 8 R e ] //3 



1 

/ r \ 2 ) 

X 

0.37 R’°* 2 
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e 

1 \6.9-10 7 / ) 



x = Downstream distance from the leading edge 

= Reynolds number = x/v 
V = Kinematic viscosity 
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A comparison of the predicted power spectrum with experimental spectra is presented 
in Figure 3. 

Cross-Power Spectra 

The final requirement in determining the characteristics of the fluctuating pressure 
field of the turbulent boundary layer is to define the ncrrow band, space correlation 
function or co-power spectral density. This parameter is the key function needed to 
describe an impinging pressure field on a structure in order to calculate the induced 
mean-square response of the structure (see, for example. Reference for the struc- 
tural response computational technique). The spatial correlation properties of a 
fluctuating pressure field can be obtained only from a careful and detailed examination 
of the field at a large number of points. Measurements by several investigators have 
shown that the co-power spectral density of turbulent boundary layer pressure fluctua- 
tions in the direction of the flow can be approximated by an exponentially damped 
cosine function, and the lateral co-spectra I density can be approximated by an 
exponential function. The general form of the cross-power spectral density is: 


S (C,1,<*>) = 4(u) A (C , >1/ <•») COS 
pp 

where A ( C, q, u) is the modulus of the cross-power spectral density, and 
0(u) is the power spectral density of the homogeneous field. 



(3) 


Here, it is assumed that the pressure field is homogeneous, in the sense that the 
cross-power spectral density is a function only of the separation distances ( £ in the 
longitudinal direction and q in the lateral direction) so that it is independent of the 
actual positions (say x and x + £ longitudinally and y and y + q laterally). Further, 
u and U are the circular frequency and convection velocity, respectively. Assum- 
ing that A ( £, q, w) is separable into its longitudinal and lateral components, and 
normalizing by the power spectral density of the homogenous field gives (Reference 1): 
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G (S,n, w) = G r (C, u))G(n, u) 
pp C n 



(£, cj) • cos 



A n (n,u) 


= C (£, u) • C (n, w) (4) 

where C (C, u) and C (n, u) are the correlation coefficients in the longitudinal 
and lateral directions, respectively. The assumed separable form leads to the pre- 
diction that the magnitude of C is constant along straight lines on the surface, 
forming a diamond pattern surrounding the origin. This characteristic is somewhat 
physically unreasonable (see Reference 2); however, for purposes of calculating 
the induced structural response the assumption of separability greatly simplifies the 
mathematics and, hence, it is generally accepted. However, Lowson (Reference 1) 
notes that a more likely form for the lines of constant amplitude would be elliptic, 
suggesting that the usual separable solution underestimates the correlation area by 
7T/2. Thus, integration of formula containing the cross-spectral density function 
should be multiplied by a factor of 7t/2 to allow for its probable underestimate of 
the correlation area at any frequency. 

Measurements of the correlation coefficients have been made by Bull and others 
(see Reference 1) and the results are presented in Figures 4a and 4b. It is seen 
that the data in Figures 4a and 4b have been collapsed based on Strouhal numbers 
C w q cj 

jj — and — — . From these data, the following empirical expressions were 

c c 

derived for the correlation coeffi Ments: 


C (C, u) = exp | -0.10 j C J u/U^ j ♦ cos j (5) 

C (q, u) = exp / -0.715 J t) j u/U c \ (6) 
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These correlation curves have been Inserted In Figures 4a and 4b for comparison. 
Typical values of the convection velocity (itself a function of frequency) fo' subsonic 
flow ere U = 0.6 U for the small scale eddies near the wall and U = 0.9 U 

C OD C 0D 

for the large scale eddies near the outer edge of the boundary layer. 


The accuracy of Equations ( 5 ) and (6 ) break down at small values of -jj- ; 

c 

however, in Reference 3 , Bull presents measured asymptotic values of the 

^ c 

correlation coefficients for small values of -g— and -g — . Based on these data, 

c r 

the Equations ( 5 ) and ( 6 ) may be corrected to include the lower frequencies, and 
the resulting expressions are: 


C ( C, u) = exp 


- 0.10 



(7) 


C (n, <*>) = exp 



• exp 



( 8 ) 


These expressions appear to be valid at both subsonic and supersonic speeds. 
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2.2 Separated Flow 


Separated flows as induced by steps, wedges, flares and other, basically rwo- 
dimensional geometric changes have undergone considerable study only in recent 
years. Considerably less dcta is available on the fluctuating p'essure environments 
within separated flow regions than is the case with attached turbulent boundary layers. 
Fur‘hermore, there are various types of separcie^ flows and little is known of the 
similarities and differences of their statistical properties. Example separated flow 
environments are listed below: 

e Blunt body-induced separation (as occur at cone-cylinder and flare -cylinder 
expc r.sion corners at subsonic Mach numbers) 

• Flare-induced, step induced, and wedge-induced separation (as occur in 
compression corners) 

• Shock-induced separation (as occur on cylinders, airfoil.., etc., beneath 
terminal shock waves at transonic speeds and due to shock wave impingement 
at supersonic speeds) 

• Boattai I -induced and rearward facing -step-induced separation (such as occur 
in the base region of launch vehicles). 

All of the foregoing environments differ to some degree in their aerodynamic structure. 
However, some basic comments can be made in regard to their fluctuating pressure 
characteristics. First, all of these environments may be reg^.ded as two-dimensional 
type separated flows having mean separation and reattachment lines which are no- u- 
to the free-stream. Second, a general characteristic is that if the flow separates 
from an expansion corner, the separation line is quite -table in that osc'llations which 
produce fluctuating pressures are not generated. However, if separation occurs, say, 
on the cylindrical portion of a payload • roud (flare induced seporrtr;n;.' . r.-f oration 
point is unstable and me/ produce significant fluctuc' : :. d pressures, particularly at 
supersonic speeds where the separation is acco:- id by an oblique shock wave. 
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Third the reottac^ment point of t*\e separated flow field produces rather large fluctua- 
ting pressure levels for virtually all type? of separated flow fields. The region within 
the separated flow field (between the separation and reattachment points) is a fairly 
homogeneous environment which is characterized by fluctuating pressure levels 
greater than those for attached flow but less than those encountered at the separation 
and reattachment points. Example data tor various separated flow fields are presented 
in the following sections. 

Overall Level 

A typical example of the fluctuating pressures resuming from blunt-body separation 
is shown in Figure 5 (results raken from Reference 18). These data were obtained 
at high subsonic Mach number or a 25-degree cone-cylinder configuration. The 

axial distribution of VpVs - uv.s a ielatively nonhomogeneous environment with 

a> 

a peak !evel which moves aft with increasing Mach number. The peak in ^ P / q 

ao 

results ..om the reattach r of of the separated flow from the shoulder. Thus, the 
extent of the separated region increases with increasing Mach number. Peak levels 
of rms fluctuating pressure reach 11 percent of free-stream dynamic pressure at a 
f r ee-stream Mach number of 0.70, and results from the instability of the reattachment 
point. It will be noted that the fluctuating pressure levels near the shoulder (X/D = 0) 
are relatively low (same order of magnitude as generally found wit. *ne homogeneous 
region of Vo-diif.ensional separated flows and typical of the environmenr for separated 
shear layerc) thus indicating that the separation point which occurs at the shoulder Is 
relatively stable* Separated flow over the boattail region of a bulbous vehicle may 
be expected to exhibit fluctuating pressure characteristics very similar to the cone- 
cylinder; however, the blunt-body separat : on on a cone-cylinder body is limited to 
the subsonic speed range, whereas, the boattail configuration may induce separation 
at all Mach numbers. 
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Typical fluctuating pressure data for flare-induced separation ere presented «n 
Figure 6 (results taken from Reference 14). These data clearly show the region of 
homogeneous separated flow, bounded on the upstream by the oscillating shock wave 

( forward peak in yj P 2 / q ), and on the downstream by the reattachment perturbations 

00 

(aft pea^ in ^ P 2 / q ). Surface fluctuating pressures for the separated flow region 
co 

range from 1 .5 to 2.7 percent of the free stream dynamic pressure. Levels associated 
with the upstream shock wave generally range from 4 to 8 percent of the free-stream 
dynamic pressure (see Reference 14); whereas, levels in the region of reattachment 
may range from 6 to 12 percent of and agree reasonably well with the reattachment 
levels for blunt body separation. Further discussion of shock-wave oscillation data 
is presented in a later section. 

The variation of fluctuating pressure level, normalized by free-stream dynamic pressure, 
with local Mach number (M^ ) for various separated flaw environments downstream 
of expansion corners is presented in Figure 7 . The regions aft of cone-cylinder 
junctures and rearward-facing steps, ond in the near woke of boattail configurations 
are represented by the data presented in Figure 7 . These environments will be 
referred to as expansion induced separated flows and it will be noted that the attendant 
fluctuating pressures exhibit the same general trend with local Mach number. The 
largest levels occured at low Mach numbers and decreased as local Mach numbers 
increased. These data represent the region of plateau static pressure and rhe tolerance 
brackets on the data represent the variations due to non -homogeneous flow within 
the region of constant static pre-.ure rather than scatter in the measurements. A good 
empirical approximation to these experimental measurements is: 

Expansion Induced Separated Flow: 



0.045 
1 + M 2 


(9) 
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This equation is similar in form to that previously proposed for attached turbulent 
boundary layers. 

Fluctuating pressure measurements for the region of plateau static pressure upstream of 

compression corners are presented in Figure 8. The regions immediately upstream of 

forward -facing ster , wedges, and flares are represented by the data presented in 

Figure 8. Also, the previous fluctuating pressure data for expansion inckiced separated 

flow, shown in Figure 7, are presented in this figure for comparison. In general, the 

compression corner data show an increase in fluctuating pressure level with increasing 

free-stream Mach number in the range, 1 .0 < M < 2.0 - reaching a constant level 

ao 

at Mach numbers above 2.0. Free-stream Mach number is used here because adeq> \ 
data is not available for determining the local Mach number in the vicinity of the 
compression induced separated flow region. Derivation of an empirical prediction 
formula for the fluctuating pressure level within compression induced separated flows 
has not been attempted at this time. 

Power Spectra 

The most comprehensive available data for power spectra of the fluctuating pressure 
within separated flows was obtained for the homogeneous region of compression 
corners c* supersonic Mach numbers (References 14, 19, 20 and 21). These data, 
presented in Figure 9, were obtained for forward facing steps, wedges and conical 
frustums. All data, represented by the cross-hatched band, showed a distinct 
similarity in spectral characteristics when comoared using normalized spectral level 
and frequency expressed as functions of local velocity, free-stream dynamic pressure, 
and local boundary layer thickne A number of velocity, length and pressure para- 
meters were used to collapse the data; however, local velocity (IJ^ ), local boundary 
layer thickness (6^ ) and free-stream dynamic pressure appeared to be adequately repre- 
sentative of the parameter dependence of the fluctuating pressures fcr the configurations 
stuaied. Power ;rru of the fluctuating pressures within the homogeneous region of 
separated flows may be represented by the following empirical formulc; 
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( 10 ) 


♦(f) u, fV<d 

/ __ <*> 
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0.83 12.15 

’T /f 0 1 

where 

f 

0 

u, 

0.17 / 

0 

L 



- 

0.045 
1 + M* 

, Figure 7, for expansion induced 
separated flows. 


■ 

the results as determined in Figure 8 for compression 
induced separated flows. 


and the subscripts l and co refer to local and free-stream conditions respectively. 

It is anticipated that Equation 10 can be used with good accuracy to predict the power 
spectra for fluctuating pr assures within the homogeneous region of expansion induced 
separated flows although it was derived based on data taken in compression corners. 

Cross-Power Spectra 

Typical cross-power spectra for the homogeneous region of two-dimensional separated 
flows are presented in Figure 10. Again, noting that the co-spectral density is the 
same as the narrow-band spatial correlation, it is seen that the separated flow exhibits 
spatial coherence very sfmila. to that of attached turbulent boundary layers. The 
damping of the sinusoidal cross spectra for separated flow is exponential at high values 
of w £ /U c as is the case for attached flow. Thus as a first approximation, the nor- 
malized longitudinal co-spectra may be represented by: 

C(C, f) = e " aU?//Uc cos iji (11) 

c 
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The damping coefficient, a, is dependent on free-stream Mach number according 

to the results presented in Reference 14. The Chyu and Hanely results show damping 

coefficients ranging from approximately 0.13 at M = 2.5 to a value of 0.33 at 

co 

M = 1.6. This suggests that the turbulence structure in separated flows decays 

OD 

somewhat more rapidly than for attached flow which has a coefficient of exponential 
decay of 0.10. It should be noted that the exponential decays represent the envelope 
of the cress-spectra for various spatial distances, ^ . Fora given value of the 
cross .pectra can be represented by the exponential envelope only at high frequencies, 
the lower limits of which increase with increasing distance between measurement 
points. 

The loss of coherence at low frequencies precludes a general collapse of the data using 
a constant damping coefficient. This problem was overcome by Coe and Rechtiei , 
Reference 20, by introducing an attenuation coefficient which is related to the nor- 
malized modulus of the cross-power spectral density by 


G 



e 


-a £ 


( 12 ) 


The normalized moduli! for available or selected transducer spacings, C , were curve- 
fitted by an exponential function using the method of least squares to obtain a non- 
dimensional attenuation-coefficient function a f ^*h in References 19-21, 


The parameter h is the height of the protuberances used to generate the separated flow 
field. Empirical approximations of the attenuation coefficient, based on the experi- 
mental results of Coe and Rechtien, are: 


a 



0.75/ in. 



< 6x 10“ 3 


(13) 
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0.75 / in. 


!i/l 

( ,S /A )o 


i- > 6x 10 -3 


l /o 


= 6 x 10" 3 


It will be noted that the longitudinal and lateral attenuation coefficients are the same 
at Strouhal numbers, f < 6 x 10~ 2 and that the lateral attenuation coefficient 

becomes larger than the longitudinal value at f 6^ ^ U ^ > 6 x 10~ 2 . It was pointed 

out in Reference 19 that this spatial characteristic indicates that the predominant 
turbulence is ncnconvective at the l.'wer frequencies and that contours of equal 


A -18 



correlation would be circular; 


whereas, at 




> 6x10" 2 the divergence 


of the longitudinal and lateral attenuation coefficients indicate a progressively 
extended correlation pattern in the direction of the free-stream with increasing 
frequency. This statement is not entirely true sine the usual separable form of the 
cross-power spectral density leads to the prediction that the magnitude of the normalized 
modulus is constant along straight lines on the surface, forming a diamond pattern 
surrounding the origin rather than a circular or elliptic pattern. Under the assumption 
of separability of the longitudinal and lateral cross-power spectra, the following 
equations (which employ the attenuation coefficient) may be used as prediction 
formula for the normalized longitudinal and lateral co-spectra. 


Longitudinal Co-Spectra 


C 



e 



$ 


cos 



(18) 


Lateral Co-Spectra 


where 



(19) 



a 



as defined in Equations 13, 14 and 15. 


= a 


as defined in Equations 16 and 17. 



2.3 


Shock-Wave Oscillation 


Generally, shock wave oscillation produces the most intense fluctuating pressure levels 
that are usually encountered by a vehicle. As for the case of separated flow, there 
are many types of shock-wave oscillation and little is known in regard to the similarities 
and differences of their statistical parameters. Typical shock waves encountered by 
vehicles are: 


• Terminal shock waves for regions of transonic flow 

• Displaced oblique shock waves as induced by the separated flow 
in compression corners at local supersonic speeds 

• Reattachment shock waves in the vicinity of the reattachment 
point for separated flows generated by both compression and 
expansion corners. 

• Impingement shock waves as caused by local bodies such as 
strop-on rockets. 

All shock waves may be expected to produce similar fluctuating pressure environments 
since the movement of the shock wave results from the interaction with the separated 

flow at the foot of the shock wave (see Reference 19) and the fluctuating pressure is 
the result of the modulation of the pressure gradient through the shock wave. A 
special case of shock wave oscillation is referred to as an alternating flow condition, 
whereby, the flow at an expansion corner intermittently fluctuates between a 
separated and attached condition. This environment is illustrated schematica'ly in 
Figure 11 fora 25 degre* me angle together with the more common terminal shock- 
wave oscillation case. Example dato for various shock wave oscillation environments 
are presented in the following sections. 
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Overall Level 


The axial distribution of rms fluctuating pressure resulting from terminal shock wave 
oscillation Is shown in Figure 11 (from Reference 18). A special case of terminal shcck 
wave oscillation results when the terminal shock wave moves forward to the expansion 
shoulder of a cone-cylinder. For this case, the flow intermittently fluctuates between 
the blunt-body separated flow condition and the attached flow condition at high sub- 
sonic (low transonic) Mach numbers. This condition represents an alternating unbalance 
between the large pressure rise through the shock wove that exceeds the values required 
to separate the flow and the small pressure rise that is too small to maintain fully 
separated conditions. 

Extremely large fluctuating pressures result from this condition; however, it should be 
noted that this phenomenon occurs over a small Mach number range and generally is of 
very low frequency. Thus for Icrge Mach number transients, this phenomenon may not 
occur. On the other hand, some experimental studies using aeroelastic wind tunnel 
models indicate that this phenomenon may become coupled with the vibrational response 
of vehicles such that flutter in the lower order bending modes would result for certain 
configurations — particularly for bulbous shaped payloads on rather slender launch 
vehicles. 

As Mach numbe* is increased above the range of alternating flow, the localized 
oscillation of the shock wave produces intense fluctuating pressures for the region in 

close proximity to the shock wave as shown in Figure 11 . The shock wave moves aft 
with diminishing strength with increasing Mach number such that the nns fluctuating 
presssure levels also decreases. In addition to the results presented in Figure 11 , 
the fluctuating pressures which occur at the separation and reattachment points for 
separated flow over compression corners (Figure 6 ) are fairly complete examples 
of shock-wave oscillation data. 
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Power Spectra 

Only recently has comprehensive data been presented on the spectral characteristics 
of shock-wave oscillation. Much of the previous data were presented in linear 
linear graphical form rather than using log— log scales. As a result, much resolution 
was lost at the high frequencies. Recent experimental data by Coe end Richtien 
(Reference 20) gives a clearly defined spectrum for shock wave oscillation at M^ = 2.0; 

however, data at other Mach numbers have not been published. Data obtained for 
three-dimensional protuberance flows do agree with the Coe and Richtien data and 
thus substantiates their limited published results. The normalized power spectra for 
shock-wave oscillation for both two- and three-dimensional protuberances (References 
18 and 19) are presented in Figure 12. The power spectrum shows a relatively steep 
roll-off starting at a Strouhal frequency (f ) of 1 x 10 2 , where the subscript 

0 denotes local velocity and boundary layer thickness upstream of the shock wave. The 
roll-off is 8 dB per octave for the range 1 x 10“ 2 < f 6^ ^ S 2 x 10 1 and above 
this range the roll-off changes suddenly to 4 dB per octave. These unique spectral 
characteristics of shock-wave induced fluctuating pressures are explained by the 
physical behavior of the shock-wave oscillation and the resulting pressure time history. 
The shock wave is basically a pressure discontinuity which becomes slightly distorted 
by the boundary layer such that a finite gradient through the shock wave is observed 
at the surface. Oscillation of the shock wave produces a wave form whl^h approaches 
a random-rectangular wave as the displacement of the oscillation increases. Superim- 
posed upon this signal is the low amplitude, high frequency disturbance associated 
with the attached boundary layer (for that portion of the signal when the shock wave 

is aft of the measurement point) and the moderate amplitude and frequency disturbances 
associated with separated flow (for that portion of the signal when the shock wave is 
forward of the measurement point). The roll-off rate of the power spectrum for a ran- 
dom-rectangular wave form is 6 dB per octave which is 2 dB lower than the experi- 
mentally observed value. Above f ^ 0 y/u o =2x 10~ ! , the power spectral density 
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for the random modulation of the shock wave diminishes below the power spectral 
density for the turbulence portion of the signal. Thus, the roll-off rate changes to a 
value roughly equal to that for separated flow since this environment is the larger 
of the two turbulence generating mechanisms (the other being attached flow). 

Noting that the power spectra for shock wave oscillation is composed of 1) low 
frequency spectral energy of the shock wave and 2) high frequency spectral energy 
of the separated flow and attached boundary layer, the resulting empirical formula 
for the power spectra may be written as a combination of power spectra of the contri- 
buting sources: 

Msw - [*< f >]l'w H + k , Ms Ma < 20 > 

where the subscripts and superscripts denote the following: 


Subscripts: 

SW 

- shock wave 


S 

- separated flow 


A 

- attached flow 

Superscripts: 

I 

- absence of viscosity (inviscid) 


H 

- homogeneous flow 


The constants, k and k are weighting functions which accounr for that portion of 
1 2 

the total energy resulting from the presence of viscous flow in the form of separated 
flow and attached flow respectively. It should be noted that the two secondary environ- 
ments (separated flow and attached flow) are not simultaneously superimposed on the 
shock wave signal but rather are time shared. This, together with the fact that these 
environments may be correla. .-d with the gross motion of the shock wave results in 

values of k and k less than 1 .0. Finally, for peak overall levels of shock wave 

1 2 

oscillation (corresponding to a point located at the mean position of the she wave) 
the contribution of attached flow is negligible in comparison to that for separated 
flow. Thus, Equation 20 may be simplified to 
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Based on the experimental data of Reference 19/ the power spectra £$(f) for 

shock wave oscillation in the absence of viscous flow normalized by local inflowing 
boundary layer thickness and velocity and free-stream dynamic pressure is given by: 



where: 




1 H 
- SW 


overall level of shock oscillation peak 
corresponding to the mean location of the 
shock wave . 



' H 
- S 


overall level of homogeneous separated flow as 
defined from Figures 7 and 8. 


The subscript 0 denotes local velocity and boundary layer thickress upstream of the 
shock wave . 
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= 1 x io -2 


'sw 


k = 0.25 

1 


determined empirically from 
experimental data of Reference 


19 


Substitution of Equations 10 and 22 into Equation 21 gives the final expression for the 
power spectra for shock wave oscillation. 
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where (fpS^/U^ *s now defined for conditions upstream of the shock wave. 

A comparison of the predicted power spectra for shock-wave oscillation with experi- 
mental measurements is presented in Figure 13. Also shown in the upper right hand 


corner of this figure is the variation in P 2 / q with distance upstream from the 
45 degree wedge. It should be noted that this prediction ferula holds true only 
at a point corresponding to the mean location of the shock wave. Oi either side 
of the shock wave, the influence of the shock dimini.' s rapidly due to its small 
displacement such that the environment is basically either attached or separated 
flow with some low frequency intermittency due to the shock wavu. 1: is convenient 
to refer to these regions as non-homogeneous c Cached at d separated flows and they 
will be discussed later in Section 2.4. 
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Cross-Power Spectra 


Very little data has been published in the form of cross-power spectra o f fluctuating 
pressures beneath oscillating shock waves. Because oscillating shook waves at a giver, 
flight condition are confined to relatively small art-s of the vehicle surface, it is. 
extremely difficult to define the spatial characteristics of the attendant fluctuating 
pressures. Fluctuating pressures in the vicinity of the shock wave are highly non- 
homogeneous; although they do appear to be related in both spectral shape and 
spatial coherence. The only significant results or. the spatial co. <»,-rence of fluctuating 
pressures in the vicinity of shock-waves are those by Coe and Rechtieri (Reference 20) 
Theii data indicate that the fluctuating pressures generated by the shock wave are 
related only at frequencies below f 6 /u = 0.08 for the region immediately 


downstream of the mean location of the shock wave (Figure 14). For the region 
immediately upstream of the shock-wave, a small degree of coh 'ence is also evident 


in this frequency range as well as ai f 6 


U >0.2. A comparison of the pcwer 
0 


spectra and coherence function shows some very interesting characteristics of shock’ 

induced fluctuating pressures.. First, the power spectra of fluctuating pressures on 

each side of the peak level point show large low frequence energy which can be 

identified as having the same basic characteristics as the shock wave spectrum for 

f 8 /U < 0.08. This is confirmed by th<* coherence of the data over the same 
0/0 

frequency range |fS ^ LJ < 0.08 j . For fS ^ U > 0.08, power speclra 

immediately upstream and downstream of the shock wave show spectral characteristics 
identical to attached turbulent boundary iayer and separated flow, respectively. Thus, 
for f6 0/ /^ > 0.8 , the spatial correlation cf fluctuating pressure immediately up- 


stream of the peak should be characteristic of attached flow; wheraas, immediately 
downstream of the peak they should be characteristic of separated flow. However, 
wh°n the *pati r rehtion is normalized by the pcwer spectral densities ro obtain 
the coherent’ funt >r, this coherence appears to be minimized due to the large 
spectrum level for the point of peak fluctuating pressure. Further discussion or. th 
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characteristic will be given later in the section on non-homogeneous attached and 
separated flows. 

The spatial deca> of the low frequency, shock induced fluctuating pressure in the 
longitudincl direction as shown in Figure 14 may be represented by an exponential 
coherence function as follows: 


y 




(24) 


A comparison of this empirical prediction with experimental data is presented in 
Figure 14. It should be noted that, as separation distance is increased, the above 
formula foils to account for the low coherence at low frequencies. However, 
because the large non-homogeneous effects associated with the flow in close proxirrvty 
of the shock wove, the application of classical statistical methods to define the spatial 
characteristics for large separation distances may be questionable. Thus* for regions 
under the peak. Equation 24 is felt to be an accurate representation of the spatial 
characteristics of the fluctuating pressures in the longitudinal direction. 


The longitudinal co-spectra may be written: 


C (C , « 



f 6 


cos 


2tt 


(25) 


Published data is not available on the transverse spatial characteristics of shock-induced 
fluctuating pressures. However, it is anticipated that these disturbances will be 
reasonably correlated over much larger distances in the transverse direction than in 
the longitudinal direction because of the continuity of the shock wave in the plane 
normal to the flow. 


C ' 2 - 
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2.4 


Non-Homogeneous A Mocked and Separated Flaws 

Non -homogeneous attached and separated flows are defined as environments which 
are basically attached or separated; however, the statistical properties of their attendant 
fluctuating pressures vary with spatial location. Examples to be considered herein are 
attached and separated flows immediately upstream and downstream of oscillating shock 
waves, respectively. The non-homogeneity may result from intermittency of the shock 
wave oscillation or from a more basic modification to the turbulence structure of attached 
and separated flow due to the motion of the shock wave. The variations in both the 
overall level and power spectra with position relative to the shock wave are evident 
in Figure 13. These data are shown in comparison with homogeneous attached and 
separated How doxa to illustrate the presence of low frequency energy due to the 
shock wave. Again, basic characteristics of the overall levels, power spectra, and 
cross-power spectra wiil be discussed for the purpose o f defining empirical prediction 
techniques for the non -homogeneous attached and separated flows. 

Overall Level 

The overall fluctuating pressure levels for attached and separated flow in close 
proximity to an oscillating shock weve are bounded on the low side by the levels of 
fluctuating pressures corresponding *o homogeneous environments and are bounded on 
the high side by the peak fluctuatiio pressures corresponding to shock wave oscillation. 

In essence, this means that the differences between the homogeneous and non-homogeneous 
fluctuating pressure levels may be attributed directly to fluctuating pressures induced by 
the oscillating shock v/ove for the case considered here. Thus, normalized fluctuating 
pressure levels for non-homogeneous flows may be defined as: 


/xF \ H /xF \ NH /xF \ H 

M A s Ur) A s (iT7 sw 


(26) 
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where the subscripts and superscripts are defined as follows: 

Subscripts A - attached flow 

S - separated flow 

SW - shock wave 

Superscripts H - homogeneous condition 

NH - non-homogeneous condition 

Under the assumption of statistical independence between the various sources, i.e., 
attached flow, separated flow and shock wave oscillation, the fluctuating pressure 
levels may be expressed as: 



where c and c are weighting functions less than 1 .0, which represent the con- 

' 2 

tribution of the shock wave to the overall fluctuating pressure level. The values of 

c and c vary with spatial location relative to the shock wave and therefore, are 
I 2 

difficult to predict. However, the above method of representation is useful in the 
prediction of power spectra for non-hor.iogeneous flows as will be shown in the next 
section . 
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Power-Spectra 

To predict the power spectra for non-homogeneous flows, a prior knowledge of the 
overall fluctuating pressure levels is required. Under the assumption of statistical 
independence between the various contributing sources, the power spectra for non- 
homogeneous environments may be written as the summation of power spectra of the 
contributing sources. Using the same symbolic representation as for the overall level, 
the power spectra for non-homogeneous environments may be written as: 





From Equations 30 and 31 , c and c are given as 

1 2 



(30) 


(31) 


(32) 
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c 


2 



(33) 


Tc determine the power spectra for non -homogeneous attached flow as caused by 
shock wave oscillation in the vicinity of the attached flow regioi.. Equations 2, 22, 
and 32 are substituted into Equation 30, which gives o form normalized by local 
conditions upstream of the shock wave: 
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Similarly, substitution of Equations 10, 22, and 33 into Equation 31, leads to the 
following expression for non -homogeneous separated flow: 



Comparison of these predictions with experimental measurements are shown in Figure 13. 
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Figure 1. Subsonic, Transonic, and Supersonic Flow Fields for Basic Vehicle Configurations 
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Figure 2. Comparison of Pressure Fluctuation Measurements beneath Attached Flows 
by Various Investigators 
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Correlation Coefficient, C (C, u) 



Figure 4a. Narrow Bond Longitudinal Space Correlation Coefficient 
for Boundary Layer Fluctuating Pressures 
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Figure 4b. Narrow Bond Lateral Space Correlation Coefficient for 
Boundary Layer Fluctuating Pressures 
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Axial Distributions of Fluctuating Pressures for Blunt Body Separated Flow; 
25 Degree Cone-Cylinder (Reference 18) 








Induced Separation (Ref. 14) 
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Figure 8. Variation of Normalized RMS Fluctuating Pressure Level with Free-Streom Mach 
Number for Expansion Induced and Compression Induced Separated Flows 
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Figure 11. Axial Distribution of Fluctuating Pressures; 25 Degree Con** Cylinder 
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Figure 12. Comparison of Power Spectra for Shock-Wave Oscillation 
Induced by Two- and Three-Dimensional Protuberances 
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O Attached turbulent boundary layer 
^ Non -homogeneous attached turbulent boundary layer 


V & Non -homogeneous separated flow 
□ Shock wove osci I lotion 

Figure 13. Longitudinal Distribution of Pressure Fluctuations and Typical Power 

Spectra in Vicinity of Supersonic Flow Separation Ahead of a 45° Wedge 
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SHOCK WAVE OSCILLATION DRIVEN BY 
TURBULENT BOUNDARY LAYER FLUCTUATIONS 


High speed aerodynamic vehicles are subject to significant fluctuating pressure environments 
from various turbulent and compressible effects. One of the more severe environments 
encountered in transonic and supersonic flow is the oscillation of shock waves, where root 
mean square pressure fluctuations are on the order of 1/10 dynamic pressure. Oscillating 
shock waves are often associated with more complex fluctuating flow fields# such as separation 
in a compression corner. Although a considerable body of experimental data enables empirical 
predictions of some fluctuating flow environments to be made with confidence, a better under- 
standing of the basic mechanisms is needed to extend available data to more general cases. 
Perhaps the most important feature of oscillating shock waves which must be determined is the 
driving mechanism of the oscillations — whether shock motion is caused by the incoming 
turbulent boundary layer, or by flowfield fluctuatings behind it. To answer this question, a 
simple model is proposed in which the incoming turbulent boundary layer drives shock motion. 
Analysis of this model shows excellent agreement with overall fluctuating pressure levels and 
spectra for shock waves in separated flow ahead of a compression corner. This analysis, 
performed under the present contract, is presented in Reference B-l and is summarized in this 
Appendix. 

Earlier treatments of the shock oscillation follow the approach of Trilling (Reference B-2) for 
interaction of a shock with a laminar boundary layer. In Trilling's analysis, frequencies were 
found for which oscillations would be self-sustaining. This would suggest that shock spectra 
would contain disorete frequencies. However, observed spectra are generally broadband, 
with no frequency peaks. The basic source of broadband fluctuations is the incoming turbulent 
boundary layer; it is therefore reasonable to look for a connection between this turbulence 
and shock motion . 

Any oscillating motion is a combination of a forcing function and a restoring mechanism. Since 

a turbulent boundary layer contains velocity fluctuations, the shock wave within the boundary 

layer will be convected upstream or downstream. A one-dimensional model of this convection 

is adopted as the forcing function, so that speed of the shock wave (relative to its mean location) 

is given by the instantaneous streamwise velocity fluctuation within the boundary layer. This 

velocity fluctuation is represented as a €p(t) , where p(t) is a stationary random function 

00 

with zero mean and < p 2 > = 1 , and € is the rms turbulent Mach number. 

The restoring mechanism is dependent on the particular flowfield geometry. Attention is 
therefore directed to a specific geometry, separated flow in a compression corner, Figure B-1 . 
The mean location of the separation point and shock wave is governed by viscous interaction 
and corner geometry. The mean location must be stable; otherwise the separation point would 
not be there. If the separation point is displaced a distance x, it will return. If x is small 
compared to separation length, it is reasonable to expect the rate of return to be proportional 
to x , so that the shock wave returns with speed -3x , where the constant 0 is a function 
of mean flow parameters and corner geometry. 
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Combining the random boundary layer convection with the linear restoring mechanism, the 
equation of motion for the foot of the shock wave is: 


u = a €|J (t) -0x 
oo 


(B-l) 


Equation (B-l ) may be integrated to give the displacement: 


= e"^ f €a f*u(£)e^d£ 
QD J 


(B-2) 


Analysis of this equation of motion leads to the following statistical properties at large time: 
• Root mean square displacement: 


<” 2 > i ■ “a,( T / 8 ) 


(B— 3) 


where = integral scale of boundary layer f'uctuations. 

The analysis leading to Equations (B-2) and (B-3) for the displacement follows 
almost exactly that for the velocity of a particle in Brownian motion with 
damping (Reference B-3). 

• Autocorrelation function: 

R w = <«(») «fr+T) > , -Br (M) 

* <x 5 > 

so that the integral scale of shock motion is T ~ 1/0. This result was obtained 
' x 

by using Equation (B-2) in the definition of R^(t). The manipulations, con- 
tained in Reference B-1 , are somewhat lengthy but straightforward. 

Fur a shock wave whose mean profile is P(x) and thickness is greater than < x 2 , the 
fluctuating pressure intensity and spectrum ere found to be: 

< p 2 > * = < x * > * (B-5) 
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4 < P 2 > 


(B-6) 



1 + 



Equation (B-5) assumes a constant gradient over the displacement distance; Equation (B-6) 
is the Fourier transform of Equation (B-4). 

The fluctuating pressure level and spectrum given by this theory were compared with experi- 
mental data obtained by Coe and Rechtien (Reference B-4). Static and fluctuating data are 
shown in Figure B-2. For the D = co , h = 8 inches case of Figure B-2 '• present theory 

predicts ( p 2 >^/q = 0.068. This is in excellent agreement with "<ta. The 

oo 

predicted spectrum is shown in Figure B-3, along with the measured spech_ >■ irom Reference 
B-4. Also shown in Figure B-3 are two spectra obtained by Robertson (Reference B-5) for 
cylindrical protuberances. 

The straightforward physical model proposed and the excellent agreement with experimental 
data lead to the conclusion that shock wave oscillation is primarily due to velocity fluctuations 
in the incoming turbulent boundory layer. The agreement of the present theory with shock 
spectra from both a two-dimensionol compression corner and three-dimensional protuberances/ 
where the incoming boundary iayer and separated shock are similar but the separated flow 
regions differ considerably, further substantiates this conclusion. 
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Figure B-l . Separation in a Compression Comer 
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Figure B-3. Comparison of Predicted Spectrum with Measured Spectra 
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MATHEMATICAL PROPERTIES OF THE EMPIRICAL SPECTRA 

Spectral data for various fluctuating pressure environments hove been analytically represented 
in the form: 


0(f) U 


00 


L 

00 


f L 
_o 

U 


(C-I ) 


1 + (f/f ) 
0 


where: U = Characteristic velocity 

L = Characteristic length 

f = Characteristic frequency 
0 

n, k = Spectrum shape factors 

The parameters f , n and k are empirically determined* Characteristic velocity and length 
0 

are chosen as those which give the most general value of f L/U over a range of flow 

0 

conditions; U is either free stream or local velocity , while L may be 8, 6*, or a body 
characteristic length. 


Equation (C-l) may be re-arranged to represent the normalized spectrum: 

i/f 


p(f> - - 7 r 

p 2 {i + (f/f) n } 


(C-2) 


/ CO — - 

0(f) df = P 2 , so that 


-oo -oo d (f/f ) 

/ *(f)df = J 

o i { 1 + ( f /n n } 
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Introduce the transformation: 


t = 
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for which d(f/f ) = — ( -j- - 1 \ n ' -L , (f/f ) 

0 n \ t / i* 0 


CD 


= r 


Aftei s light - manipulation. Equation (C-3) becomes 

.1 


-co , -i l-l k- 1-1 

/ V(f) df = -L / (l-t) n t dt 

0 0 


= 1 


The definition of the Beta function is (Equation 6.370 o? Reference C-»): 


B(x,y) = f t X '' (1 - t) y dt. 


Re (x) > 0 
Re (y) > 0 


so that Equation (C-4) is: 



n 




= 1 


The Beta function may be represented in terms of the Gamma function (Equation 6.381 
Reference C-l): 


B(x,y) 


r(x) r(y) 

1 W 


so that the present condition is: 
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Noting the well known properties of the Gamma function; 


ro) = i 

r(x+i)= xr(x) 
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(C-8) 


(C-9c) 

(C-9b) 



Equation (C-8) Is satisfied by: 


k = 1 + — (C-10) 

n 


It is obvious from Equation (C-3) that for each value of n there is a unique correct value of k. 
Noting that Equation (C-10) and Its inverse are single valued. Equation (C-10) is therefore the 
unique solution to Equation (C-3). 

With Equation (C-10), the three free parameters f , n and k are reduced to two: f and 

0 0 

either n or k. The prediction formulae presented in the body of this report and in Appendix A 
did not make use of this relation, but rather fitted all three parameters independently, using 
three points in each measured spectrum. The values of n and k obtained are listed in Table C-l . 


TABLE C-l 

VALUES OF n AND k USED IN EMPIRICAL SPECTRA 


Environment 

n 

k 

Attached Flow 

0.9 

2.0 

Separated Flow 

0.83 

2.15 

Shock Oscillation 

1.55 

1.7 


Figure C-l shows these values of n and k as compared to Equation (C-10). The agreement 
is quite good. Although it should be understood that the present point is strictly a mathematical 
manipulation, the agreement shown in Figure C-l provides additional justification for the use 
of the form of Equation (C-l ). 

The values of n and k shown in Table C-l are sufficiently close to agreement with Equation 

(C-10) that they are retained unchanged in the prediction formulae. If at some future time it 

is desired to revise n and k, this can be easily done. It should be noted that the most natural 

way of fitting the two independent parameters, f and n or k, is to utilize data at low 

0 

frequencies ( f <3C f ) to obtain f Q , and high frequency data (f»f Q ) to obtain n and k. 

At f » f Q , the slope of the spectrum on a log-log plot is - nk . 
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Figure C-1 . Comparison of Fitted Values of n and k with 
Equation (C-10) 
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WYLE LABORATORIES 


COMPUTER PROGRAM DESCRIPTION 


Program Number: 
Authors: 

Date : 

Source Language: 
Computer: 


71/004S-1 A 
D .Lister, K .Plotkin 
May 1973 
Fortran IV-H 
XDS Sigma 5 


1 .0 PROGRAM TITLE 
Power Spectral Density 

2.0 PURPOSE 

Identical to Program 71/004S-1, written by D. Lister, except that the present version 
is a subroutine and does not provide plots. 

3.0 METHOD 

Same as Program 71/004S-1 . 

4.0 COMPUTER CONFIGURATION 

The required hardware is: XDS Sigma 5 Computer with 16 K core, card reader and 
line printer. 

5.0 CALL SEQUENCE 

The program is entered by the call: 

CALL PSDIN (FB, FE, ANN, TIT, IABCD, QS, DB, UL, PSQUAR, 

FZERO) 

All input variables correspond to those described for 71 /004S-1 , with the following 
additional notations: 

• IABCD is a 4 element array, such that 

IABCD (1) = IA 
IABCD (2) = IB 
IABCD (3) = IC 
IABCD (4) = ID 
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• PSQUAR is a 3 element array, such that 

PSQUAR (1) = PQSA 
PSQUAR (2) = PQSB 
PSQUAR (3) = PQST 

0 FZERO is a 3 element array, such that 

FZERO (1) = FDUA 
FZERO (2) = FDUB 
FZERO (3) = FDUT 

After a call to PSDIN has been made, subsequent entries for identical values of FB, 
FE and ANN may be made by the statement 

CALL PSD (TIT, IABCD, QS, DB, UL, PSQUAR, FZERO) 

6.0 OUTPUT 

The output is identical to 71 /004S-I . 
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WYLE LABORATORIES 


COMPUTER PROGRAM DESCRIPTION 


Program Number: 
Author: 

Date: 

Source Language: 


73/004P - 1 
K . J . Plotkin 
April 1973 
Fortran IV - H 


1 .0 PROGRAM TITLE 
Correlation Function 

2.0 PURPOSE 

Given the local velocity and boundary layer thickness, the program computes and 
prints narrow band cross correlation functions for attached turbulent boundary layer 
flow, separated flow, and protuberance induced separated flow. 

3.0 METHOD 

The prediction method of Reference 1 is used. The tormulce, presented in Reference 
1 , for longitudinal and lateral correlation functions may be written 


Longitudinal: 


f) - exp(-|«l/L.) cos 2*%- 
’ c 


Lateral: 


A yj = exp (- i^l/L^ ) 


where £ and rt are streamwise and transverse coordinates, f is frequency and u 

c 

is narrow band convection velocity. The present notation differs slightly from that 
of Reference 1 in that A^ here includes the cosine term . The notation L^ and 

L^ is introduced here, with these two quantities representing the exponential decay 

length scales of the correlation functions. 


The quantities L ^ 


L -q and 


u are functions of frequency. 


The program computes L^ , L^ , A^ and A^ for attached turbulent boundary layers 

and flare induced separated Hows. Longitudinal quantities. La and A* , are com- 
puted for protuberence induced separated flow. 
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4.0 COMPUTER CONFIGURATION 

The required hardware is: XDS Sigma 5 Computer with 16K core, card reader, and 
line printer . 

5.0 CALL SEQUENCES 

This program is written as a subroutine, with four entry points. The entry names and 
call sequences are as follows: 

Initial Entry: Initializes frequency range and computes decay length scales. 


CALL CLENFR (FB, FE, ANN, DB, UL, IATT, ISEP, IPROT, IGNORE) 


where 


FB 

= 

Lower limit of frequency range of interest 

FE 

= 

Upper limit of frequency range of interest 

ANN 

= 

Number of frequency points per decade 

DB 


Local boundary layer thickness, feet 

UL 

= 

Local velocity, feet second 

IATT 


j I if computation desired for attacned flow 
l 0 if not desired 

ISEP 

- 

{ 1 if computation desired for separated flow 
lo if not desired 

IPROT 

- 

11 if computation desired for protuberance flow 
(0 if not desired 

IGNORE 

:= 

-1 


An entry to CLENFR compares L^ and L^ for the desired frequency range and 

saves these for later printout. Storage is in a COMMON block described in Sec- 
tion 6.0 . Any or all of the environments mcy be selected for computation . The 
first call to the routine, and any call with a change in FB, FE or ANN, must be to 
CLENFR. Each time CLENFR is called, the frequency range is printed out. 


Subsequent Computation of L. and 


If it is desired to compute L^ and 


’»7 ' 


L^ for different flow conditions over the same 
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frequency range as before, this is accomplished by: 


CALL CORLEN (DB, UL, IATT, ISEP, IP*^T) 


Computation of A^ and 


A 


r) : 


A^ and A^ are computed as functions of distance for a specified frequency. The 
call is 


CALL CORFUN (DB, UL, IATT, ISEP, IPROT, IGNORE, 
IFREQ, XB, XE, DEL) 


where 


IGNORE 

IFREQ 

XB 

XE 

DEL 


1 if previously computed values of L. and L^ 
are to be used 

C if L^ and Ljj are to be re-computed for the 
speci ‘ied values of DB and UL 


Index of the frequency for which 
to be computed 


A^ and 


are 


Lower limit of separation distance range of interest 
Upper limit of separation distance range of interest 
Increment between distance points to be computed 


Longitudinal and lateral correlation function are computed for the specified separa- 
tion distance range and the frequency corresponding to IFREQ. These values are 
saved for later printout. 


Printout of Computed Values: 


CALL CORPRT (TIT, ILEN, ICOR) 


where 


TIT = Name of a 20 element alphameric array containing 
a title to be printed out 
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ILEN 


I COR = 


'1 if l* and L« are to be printed 
\0 if not 

( 1 if A. and A„ are to be printed 
(0 if nor 


A call to CORPRT with ILEN = 1 results in a printout of the lest computed values of 
L^ and L^ for all three environments as a function of frequency. Strouhal number, 

u. , is also printed . 

L 

A call to CORPRT with ICOR = i results in a printout of the last computed values of 
A^ and A^ for all three environments. 

The title TIT is printed above each table. Also printed are parameters of the fre- 
quency range . 

A sample output is cttached . The printout is sufficiently annotated so as to be self- 
explanatory. 

6.0 STORAGE AND RETRIEVAL OF COMPUTED RESULTS 

In the event that access is required to values computed by this suLrout'ne , computed 
values are stored in a COMMON block designated COLD. To gain access, the 
following statement must be inserted into the main program or other subroutine: 



COMMON /COLD/ F, ST, ATTLX, ATTLY, SEPLX, SEPLY, PROTLX, 
LTATT, UCSEP, X, CORFAX, CORFAY, CORFSX, SORFSY, CORFPX 


where each name is that of a 100 element one dimensional array. Each array con- 
tains the foliowing*. 


F 

ST 

ATTLX 

ATTLY 

SEPLX 

SEPLY 

PROTLX 

UCATT 


Frequency 

Strouhal number, f<5/i^ 

L^ for tached flow 

L^ for attached flow 

L^ for separated flow 

L^ for separated flow 

L for protuberance ffow 
x 

Narrow band convection velocity for attached flow 
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UC5EP 

X 

CORFAX 

CORFAY 

CORFSX 

CORFSY 

CORFPX 


Narrow band convection velocity for separated and 
protuberance flow 

Distance coordinate for correlation function 
A^ for attached flow 
for attached flow 
A^ for separated flow 
A^ for separated flow 
A * for protuberance flow 


Different names may be used in the COMMON statement; however, the order must 
be preserved exactly and each array must be dimensioned exactly 100. The name 
COLD may not be used for any other purpose in the program . 


Because of the dimension of 100, the frequency and distance parameters must be such 
that r.o more than iOO points are computed for each . If the frequency parameter 
(FB, FE, ANN) do not satisfy this restriction, CLENFR will substitute an acceptable 
value of ANN and proceed . If the distance parameters (XB, XE, DEL) do not satisfy 
this restriction, a message will be printed and the job terminated. 


7.0 REFERENCES 

1 . Robertson, J. E., "Prediction of In-Flight Fluctuating Pressure Environments 
Including Protuberance Induced Flow," Wyle Laboratories Rerearch Staff 
Report WR 71-10, March 1971. 
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WYLE LABORATORIES 


COMPUTER PROGRAM DESCRIPTION 


Program Number: 
Author: 

Date: 

Source Language: 


7 3/005P-1 
K. Plotkin 
May 1973 
Fortran IV-H 


1 .0 PROGRAM TITLE 
Transonic Flow 

2.0 PURPOSE 

To compute fluctuating pressure environments on the cylindrical portion of c cone 
cylinder at transonic Mach numbers 0.6 £ M <0.98. 

3.0 METHOD 

The program uses the prediction schemes presented in Reference 1 thru 3 for overall 

fluctuating pressure levels and terminal shock wave oscillation. There are three 

basic cases, as a function of Mach number and cone angle. For M < .6, attached 

flow fluctuating pressure level is computed from Equation 1 of Reference 1 . For 

6<M<M (M = attachment Mach number), overall level is computed from Steven's 
o a 

prediction scheme (Reference 2) as cited in Reference 3, For attached flow and 
shock wave oscillation, overall level and shock location (relative to cylinder dia- 
meter) are computed from data presented in Figures 13 thru 16 of Reference 3. 

4.0 COMPUTER CONFIGURATION 

The required hardware is: XDS Sigma 5 Computer with 16 K core. 

5.0 CALL SEQUENCE 

The program is entered by a call : 

CALL TRANS (THETA, EMINF, XOVERD, CP, NTYPE) 

where 

THETA = Cone half angle 
EMI NF = Free stream Mach number 
XOVERD - Shock wave location divided by cylinder diameter 
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The program is entered with values of THETA and EMINF. Values of XOVERD, CP 
and NTYPE are returned. NTYPE may attain the following values: 


0 : M < 0 .6 attached flow exists and has been 

ao 

computed. XOVERD is set equal to 0. 

1 : 0.6 < M < M . Separated flow exists and has 

ao a 

been computed . XOVERD is set equal to 0. 

2 : M < M < 0.98. Attached flow and shock 

a ao 

oscillation exists and has been computed. 

3 : An unacceptably high value of M (>0.98) 

has been specified. 


6.0 INPUT AND OUTPUT 

This subroutine contains no READ or WRITE statements; values returned are to be 

hcndled by a main program provided by the user. 

7.0 REFERENCES 

1. Robertson, J. E., "Prediction of In-Flight Fluctuating Assure Environments 
Including FVotuberance Induced Flow," Wyle Laboratories Research Staff 
Report WR 71-10 , March 1971. 

2. Himelblow, Harry, "Aeroacoustic, Vibration and Shock Environments for the 
Space Shuttle Orbiter," Paper presented at Space Shuttle Dynamics and 
Aeroelasticity Technology Working Group Meeting, Ames Research Center, 
Moffet Field, California, November 8-9, 1971. 

3. Robertson, J. E., "Preliminary Estimates of Space Shuttle Fluctuating Pressure 
Environments," Wyle Laboratories Research Staff Report WR 72-10, August 
1972. 
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WYLE LABORATORIES 


COMPUTER PROGRAM DESCRIPTION 


Program Number: 
Author: 

Date: 

Source Language: 
Computer: 


71/ 004S-1 
D.M. Lister 
July 1971 
FORTRAN IV-H 
XDS Sigma 5 


1.0 PROGRAM TITLE 
Power Spectral Density. 

2.0 PURPOSE 

Given the local velocity, the free-stream dynamic pressure, the local boundary 
layer thickness and the frequency range of interest, the program prints and plots, 
the pressure spectral density, the normalized pressure spectral density and the 
fluctuatirq pressure level. 

3.0 METHOD 

The methxl employed is given >n Section 3. 2. 2. 4, page 31 of Reference 1 . 

4.0 COMPUTER CONFIGURATION 

The required hardware is: XDS Sigma 5 computer with 16 K core, card reader, 
line printer and calcomp plotter. 

5.0 PUNCHED CARD INPUT 

This is best illustrated in tabular form as shown on the following page. 
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Variable 

1 

1 

Description 

Format 

Columrn 

f b 

FB 

Initial frequency of interest 

F10.0 

I -1 0 

f e 

FE 

Final frequency of interest 

F10.0 

11-20 

npd 

ANN 

Number of frequency pt^decade 

F10.0 

21-30 


XDECADE 

Length in inches of decade of n axis 

F10.0 

31-40 


YDECADE 

Length in inches of decade of y axis 

FIO.O 

41-50 


YTDB 

Length in inches of 10 dB on y axis 

F10.0 

51-60 

Title 

n(i), 

1 = 1, 20 

A Title cord 

20A4 

1-80 




I 1 

5 


IB 

1 

I I 

10 


1C 

/ Zone parameters 

I 1 

15 


ID 

) 

1 I 

20 


IP1 

) 

I 1 

25 


IP2 

Plot parameters 

I 1 

30 


1P3 

I 

I 1 

35 

q ao 

QS 

rree-stream dynamic pressure 

FIO.O 

1-10 

5 / 

DB 

Local boundary layer thickness 

FIO.O 

11-20 

U 


Local velocity 

FIO.O 

21-30 


. QSA 

Normalized mean square FPL 
Attached Flow 

FIO.O 

1-10 

W U l\ A 

FDUA 

Normalized characteristic frequency 

FIO.O 

11-20 


Attached Flow 



(W) R 

PQ SB 

Normalized mean square FPL 

FIO.O 

1-10 

0 0 o 


Separated Flow 



%V"t\ j 

FDUB 

Normalized characteristic frequency 

FIO.O 

11-20 


Separated Flow 



(P*/q*L 

POST 

Normalized total mean square FPL 

FIO.O 

1-10 

00 1 


Shock Oscillation 



(! „V U / >c 

FDUT 

Normalized characteristic frequency 

FIO.O 

11-20 


Shock Oscillation 




PQSB 

1 Normalized mean square FPL 
Separated Flow 

FIO.O 

21-30 
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Note 

Card 

# 

Variable 

Fortran 

Symbol 

Description 

Format 

Columr ; 

9 

8 


POST 

Noimalized total mean square FPL 
Shock Oscillation 

FI 0.0 

1-10 

9 

8 

(f .vv c 

FDUT 

Normalized characteristic frequency 
Shock Oscillation 

F10.0 

11-20 

9 

8 

1 

PQSA 

l 

Normalized mean square FPL 
Attached Flow 

F10.0 

21 -JO 

I 


NOTES: 

1 


The various output quantities are computed for values of frequency lying between f, 

th ® 

and • The I frequency point from Is given by the formula: 


log (f. ) = log (fj + ~ log (10) 


2 For third octave band frequencies, Input npd = 10, 

3 XDECADE Is used by all three piots. It Is suggested that XDECAOE 1$ selected such 
that XDEC aDE multiplied by the number of x decades Is less than ten. 

(e.g., for frequency ranr,e 1 Hz to 10 4 Hz then XDECADE Is best chosen to be two.) 

4 The program inputs this card and prints Its contents prior to the main output table, 
(i.e., a title card.) 

5 At the end of processing the program returns to rend a new title card and a new 
parameter card. If the sum of the zone parameters is zero a new card type one 1$ 
read. If FE is not greater than FB then the program stops. 

6 The zone parameters must be equa‘ to zero or one. A one indicates that the structure 
is subjected to the indicated flow condition (see Reference 1, page 32). Possible 
combinations of the four zone parameters are: (1 ,0,0,0), (1 ,0,0,1 ), (0,1 ,0,0), 
(0,1, 1,0), (0,0, 1,0) and (0,0, 0,1). 

7 The plot parameters must be zero or unity. A unity indicates that a particular 
plot is required. IP1 = 1 then a plot of normalized pressure spectral density against 
frequency (log/log) is obtained. IP2 - 1 then a plot of pressure spectral density 
against f r equency (log/log) is obtained. IP3 - 1 then a plot of fluctuating pressure 
level (dB; re: 0.0002 dynes/cm*) against frequency (lin./log) is obtained. 
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8 


Units of these variables most be consistent. Note FPL is referenced assuming 
ft. lb. sec units and P ^ equal 0.0002 dynes/cm 2 , 

9 If I A = 0 then no card type five required. 

If IB =0 then no card type six required. 

If 1C = 0 then no card type seven required. 

If ID = 0 then no card type eight required. 

6.0 OUTPUT 

The line printer output is adequately annotated. 

7.0 STORAGE 

Just over a quarter of the available core is used. 

8.0 REFERENCE 

Robertson, J.E., "Prediction of In-Flight Fluctuating Pressure Environments Including 
Protuberance Induced Flow," NASA CR 119,947, March 1971 (N71-36677). 
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TYPICAL COMPUTER PRINT-OUY 
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TYPICAL COMPUTER PRINT-OUT (continued) 
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VVVLE LABORATORIES 

computer RR06RAM 0T5£*iption 


Program Number: 
Author: 

Date: 

Source Language: 
Computer: 


73/003P-1 
K. J. Plotkin 
February 1973 
FORTRAN 1V-H 
XDS Sigma V 


1 .0 PROGRAM TITLE 
Plume Separation 

2.0 PURPOSE 

Given axisymmetric launch vehicle and trajectory parameters, the program computes 
properties relative to plume-induced separation. This includes separation length and 
plateau pressure, plus various properties of the plume and parameters describing the 
coupling to the unsteady pressure field. 

3.0 METHOD 

The analytic models incorporated in this program ore described in detail in vVyle 
Laboratories Research Staff Report WR 73-3. 

4.0 COMPUTER CONFIGURATION 

The required hardware is: XDS Sigma V computer with 16 k core, cord reader, and 
line printer. 

5.0 PUNCHED CARD INPUT 

Vehicle and trajectory parameters are input in metric units. Specific format for each 
card is as follows: 


Card 1 


Descriptive title of vehicle. Up to 80 alphameric characters 


Card 2 


Columns 1-10 
Columns 1 1 -20 
Columns 21 -30 
Columns 31 -40 
Columns 41 -50 
Columns 51 -60 


Vehicle radius, meters 

Total thrust Newtons 

Plume drag to thrust ratio 

Exhaust gamma 

Exhaust Mach number 

Engine chamber pressure, N/m 2 


Format for each item on this cord is FI 0.0. 
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Card 3 — Descriptive title of trajectory. Up to 80 alphameric characters 

Card 4 — N, H » number of trajectory points to be read and computations 
performed for. Format 1 5. 

Cards 5 through 4+N — Trajectory parameters 

Columns 1-15: Time After Lift-Off, seconds 
Columns 16-30: Altitude, kilometers 
Columns 31-45: Flight Mach number 
Columns 46-60: Free-stream dynamic pressure, N/m 2 

Format for each item on these cards is E 15.10. 


6.0 OUTPUT 

The first items output are a table of atmospheric properties and a table of separation 
conditions contained in rhe program. This is followed by output of the vehicle 
parameters contained Card 2, prefaced by the title input in Card 1 . The title 
input in Card 3 is then printed, followed by a listing of the trajectory parameters 
input. Five (5) tables of computed properties are then printed. Each of these 
tables consists of the time after lift-off, followed by various properties. The 
printout is sufficiently annotated so as to be self-explanatory; the column headings 
are FORTRAN representations of various expressions defined in WR 73-3. A sample 
output list is presented herein. Written below each column is the exact notation 
used in WR 73-3, or a definition for those items not specifically denoted in the report. 

It should be noted that the program does not explicitly indicate whether or not the flow 
is separated. The output listing must be examined to determine this. The computed 
values represent separated flow when the plume angle is greeter than the separation 
angle. In the sample output, flow is not separated at 84 seconds, but is at all other 
computation times. 


D-24 



ATf'OSPKfc.wt TAbLF 


cccccccc 

♦ ♦♦♦♦♦♦♦ 

o 0 in 0 *-*000 

ON « CO 00 O 4- 4 

a c<ro'C' 0 >ox(r 1 
u C0000C00 
\r. c 0 0 0 c o 0 <r. 

c<rcr(T^co\ 0 ' 
i o<u (u o 4 4 
O 'Si 4 ^ <o o m n 

CCCCCCCC 




w-t 

o 

o 

O c 

O 

o 

c 

O 

*1 

»-• r* *4 


ru 

ft! 

•X' 

ft; 

ft 

ft! 

fti 

fti 

fti 

0 

0 

0 

0 

0 


x-f 

^4 

•»* 




fti 


c 

c 

c 

c 

o o 

O 

o 

c 

o 

c 

c o c c c 

c 

O 

c 

o 

c 

c 

o 

o 

o 

o 

c 

O 

o 

c 

c 


c o 

o 

o 

o 

o 

c 

c 


♦ 

♦ 

♦ 

♦ 

♦ ♦ 

♦ 

♦ 

♦ 

♦ 

• 

1 1 1 1 I 

• 

• 

• 

t 

• 

t 

• 

ft 

• 

• 

• 

• 

1 

• 

• 


♦ ♦ 

♦ 

♦ 

4 

4 

4 

4 


U‘ 

U* 

ui 

U* 

Ui Ui 

Ui 

Ui 

Ui 

ui 

Ui 

U Ui U lii Oj 

ui 

U 1 

UI 

Ui 

ui 

Ui 

Ui 

Ui 

u; 

Ui 

ui 

Ui 

Ui 

Ui 

Ui 


ui ui 

Ui 

ui 

ui 

Ui 

UI 

Ul 


00 

00 

.T 

rs 

vC IT. 

IT 

■4 

*4 

W -4 

x 

30 4 4 (V Tv 

X 

0 


0 J 

IT 

0 

X s 

X. 

c 

0 

X 

0 

4 

>0 

X 


c c 

0 

>0 

cj 

X 

0 

X 


-4 

-4 

4 

0 i 

fs X 

IT 

|s 

0 

OC 

or 

fs 0 C 0 C 0 C IS 

X 

ft. 

4 

r. 

0 

0 

00 

00 

0 

ps 

ps 

ac 

IT 

X 

c 


c c 

X 

*- 

c 

c 

0 

0 


IT 

0 

0 

oi 

O' <T 

0 . 

0 

0 

0 

00 

oO 0 0 0 0 

0 

rs 

rs 

P's 

rv 

0 

0 

0 

0 

0 

ps 

oO 

X 

is 

0 

a 

c c 

0 

c 

c 

0 

0 

0 . 

u. 

(Tv 

O' 

0 

01 

01 O' 

01 

01 

01 

01 

0 

0 0 0 0 0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

0 

0 

0 

0 

0 

0 

Ui 

o o 

0 

0 

o 

0 

0 

0 


01 

u* 

-4 

o 

ao o 

0 ) 

XI 

x 

4 

0 

or\ is x x> 0 

4 

n 

J0 

X'* 

x 

rs 

X 

4 

w 

X 

Ui 

X 

«-i 

ai 

*-» 


CJ o 

0 

0 

o 

0 

0 

4 ) 

•— 

-4 

< 

(T 


rs ir 

0 

X 

4 

* 

c 

C^lTNC 

4 

0 * 

is 

X 

IT 

4 

X 

w~ 

X 

X 

0 

0 

X 

X 

0 


c c 


0 

c 

0 

X 

X 

© 

ft 

o 

0 

o 

X) 0 


Pn. 

X 


0 

X* 0 ft. C 4 

C 

X 

if) 

X 

0 

0 

X) 

fs 

•*- 

ft* 

O 

«i 

rs 

c 

in 

0. 

o c 

*-• 

PS. 

c 

O 

0 

X 

X 

ft 

c 

«-• 

X 

ft «- 

*- 

ft* 

x 

<V 

OC 

4 X 4 IT X 

w- 

X 

a 

0 

0 

0 

ft. 

t— 

r 

c 

C 

0 

0 

0 

c 


c c 

ft 

X 

c 

|S 

ft 

*» 

_r 

H 

— 

flC 

vC 

X -4 

0i 

ft. 

^4 


0C 

X 4 0 (V<- 


0 

rs 

.0 

0 

ft* 

fti 

— 

w* 

— 

X 

X 

4 

0 

0 


fti ft 

ft. 

ft 

-r 

4) 

X 



c 

o 

o 

o 

C Q 

O 

o 

Q 

o 

o 

O C C C O 

o 

c 

c 

# 

O 

c 

o 

c 

o 

c 

o 

o 

O 

c 

Q 


3 a 

o 

o 

o 

c 

c 

o 


Ui 

© 



n 

•r 

if) 

41 

n 

.ft 

IT 

:ft 

4) 

4 

4 

4 

4 

4 

4 

4 

0 

m 

0 

rr 

0 

0 

0 

0 

0 

ft 

ft 

ft 

ft 

ft 

ft 

U 


0 


rs 

c 

o 

c a 


o 

c 

o 

a 

c 

C 

O 

C 

O 

o 

O 

O 

c 

C 

O 

O 

c 

C 

C 

C 

c 

c 

O 

O 

C 

o 

c 

O 

O 

C 

o 



C C 

c 

C 

c 

3 

c c 


4 

4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

♦ 

4 

4 

4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

♦ 

• 


4 4 

♦ 

4 

4 

4 

♦ 4 


al 

Ui 

Ui 

ul 

Ul 

Ul 

til 

ui 

Ul 

Ul 

Ul 

Ul 

Ul 

Ul 

Ul 

UJ 

Ul 

Ui 

Id 

ui 

UJ 

41 

UJ 

Ul 

Ul 

ui 

ui 

Ul 

Ui 

Ul 

ui 



UJ Ul 

ul 

ul 

Ui 

ui 

ui ui 


Q 

3 

ft 

a 

42 

ft 

-J 

4 

0 

4 

0 

0) 

n 

41 

0 

X 

X 

0 

0 

f) 

0 

X) 

s. 

XI 

X) 

ft 

0 

ft 

X 

4) 

£ 

v0 


X) X 

X 

*-« 

X 


4 1 4 


O 

to 

o 

PS 

«u 

4 

in 

j0 

rs 

0 

4) 

in 

ft 

4 

O 

0 

X 

X 

CM 

0 

0 

«-• 

4 

X 

X 

ft 

rv 

0 

4) 

X 

•-* 

z 

a. 

4)0 


0 

0 

4) 

X 4) 


o 

rs 


4) 

O 

ft 

rs 

0 

>. 

c 

4 

N 

c 

4> 

4 

ft 

0 

4) 

r* 

X 

fs. 

4 

4 

4 

X 

ft 


0 

ft 

0 

C 

o 

u> 

0 0 

0 

0 

0 

0 

0 0 

u. 

o 

0 

4 

4 

X 

ac 

0 

0 

rv 

4 -* 

O 

OO 

4 

4*. 

0 

0 

ft 

O 


fs 

0 

4) 

41 

o 

X) 

o 

0 

'C 

ft 

ft 

4) 


in 

'n 0 

0 

0 

0 

0 

0 0 

*? 

c 

^4 

o 

X 

*-» 

0 

0 

o 

ft 

ft 

0 

♦ 

r> 

0 

»-• 

b 

X 

W+ 

ft 

0 

4 

41 

0 

4 

0 

0 

X 

*-* 


1^ 

C 

►- 


Ci 0 

0 

0 

0 

0 

0 0 

•— 

o 

o 

x 


ll) 

0 

0 

rs 

ifi 

0 

0 

fs* 

*- 

X 

4* 

rs 

o 

4 

IT 



0 

4 

0 

ft 

PS 

ft 

0 

41 

0 

42 


X 

0 0 

0 

0 

0 

0 

0 0 

a 

a 

in 

x 

ft 

X 

4 

0 

*-• 

0 

4) 

ft 

4 

fS. 

X 

r* 

0 

0 

0 

X 

rs 

rs 

0 

0 


ft 

fs 

ft 

4 

4) 

w* 

4 

o 


— 4 

0 

X 

X 

ft 

* 0 


fti 

0 

«-» 

ps 

in 

X 

0 

4 

o 

4) 

in 

o 

0 

f-4 

4) 


30 

41 

0 

N 

00 

«-« 

41 

0 

o 

0 

ft 

X 

PS. 

0 

ft 

z 


41 ft 

0 

C3 

ft 

X) 

PS * 


x 

N 

X 

4 

0 

ft 

w-* 



rs 

in 

4 

ft 

ft 



X 

X 

4 

0 

ft 

ft 

•-* 



Pv 

4? 

4 

0 

ft 

ft 

© 

i 

ft ft 


ft 

ft 

ft 

ft 0 


• 

• 

• 

• 

• 

• 

• 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

i- 

• * 

• 

* 

4 

• 

• • 


coooooooooooooouoooooooooooooao 



O J> 

o 

O 

o 

o 

O O 


z 

© 


• • 

o o 

• 

o 

• 

o 

• 

O 

• 

o 

4 

o 

• 

o 

4 

o 

4 

o 

4 

o 

4 

<J 

4 

o 

4 

CJ 

4 

o 

4 

o 

4 

o 

4 

o 

4 

o 

4 

o 

4 

o 

4 

CJ 

4 

CJ 

4 

o 

4 

CJ 

4 

o 

• 

o 

4 

o 

4 

o 

4 

o 

4 

o 

J 

CJ 

a 

o 

o 

o 

o 

o 

o 

o 

o 

CJ 

CJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CJ 

o 

o 

o 

o 

o 

o 

o 

CJ 

o 

CJ 

CJ 

CJ 

o 

o 

u 

CJ 

CJ 

o 

CJ 

CJ 

o 

CJ 

o 

o 

CJ 

CJ 

CJ 

o 

CJ 

CJ 

CJ 

CJ 

o 

CJ 

CJ 

CJ 

ft 

4“ 

41 

X 

o 

ft 

4 

4) 

X 

o 

ft 

4 

4) 

X 

CJ 

ft 

4 

X 

X 

CJ 

ft 

4 

X 

X 

o 

ft 

4 

X 

X 

o 






^4 


44 

44 

ft 

ft 

ft 

ft 

ft 

0 

0 

0 

0 

0 

4 

4 

4 

4 

4 

m 

«n 

4> 

tf) 

in 

X 


-J 

< 



h(u n 4 i) <o rv w 


D-25 




o 

n o 

o o 

m o 

O f*> 

• 

in • 



D-36 



0»'5b999994F*n? (1» 20999994 3E+C1 r,» 358fU0nCCE +C5 

0»212999M/8t+0? n»27S99?9?8E+(j' C « P427 7OOCCE+05 


IT -# * * ■* -* 

tz u c: o o o 

♦ ♦♦♦♦♦ 

U W ut Ul lit 

cocccc 
o u o o o o 
c c c c c c 
c c c o c c 
c c c c c c 

4 j f) (Tv O 0 \ CM 
AHtmOMO 
* A. N 4* o m 
*- cr ir <r fv *- 
»»•••• 
t c c o u c 


o e o o o o 

♦ ♦♦♦♦♦ 
ui ui ui ui ui ui 
m r»> a o vO o or* 
n ^ c ^ o ^ 
tr ^ (r (r c O' 

O' 0^ 0> (T C 0> 

O' O' O' o' o <r 

O' O' O' O' O O' 

«v O' c r* tc 

iTCNDOJC 
(*> * in *+ 


c c c c c c 


A. (V (V. Al (VI (1 

c c o o c c 

♦ ♦♦♦♦♦ 

ll J uJ ui ul Ui 
O OOO'O' O' 
o o o m >o 

COCOO'O' 
CCCO'O'O' 
rj y o o' o' O' 
O O J ON ON O' 
O O ii W flO ^ 
co m jc ♦ o « 

cw <n *n *■ m w 

• ••••• 

o c c c o o 

< 

£ 


J •••••••• 

< ^OuOvO^AJOO 

Z ao O' O *-• (Vi r*> ■#■ -*■ 

• ij «H +“% H Ti H 

jc r 

© 

z •* 

© 

-I 

-J 

© 

£L 

* 


D-27 



Tim. PLATE AU PrttSSUKf. SFE'AKAT IBN ANKLET PLUPF ANKLF PLUME RADIUS SFPAKATUSK LFNHTP 


— o 

*-• 

oj 

CVI 

(Vi 

(V 

CVI 

CJ u 

u 

u 

u 

O 

c 

u 

4 4 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

U UJ 

u 

u. 

Ui 

lx 

u. 

U 

>om 

rr* 


4* 

>0 

rv 


f*“ 

c; 

rv 

IT 

<r 

4 

* 

^ 4 

r* 

c 

ao 

rr 

<r 

•-* 

•“»(U 

u 

oi 

rv 


0J 

00 

c X 

><; 

r- 

Cr 

c 

aC 

4 

4 ev- 

Ot 

rv 

rv 

c *i 

n 

4 

er 4 

ac 

vC 

4 

in 

w* 

00 

x tr 

rr 

u 

rv 

in 

rr 

(C 

1*) 

in 


»— 

Oj 

OJ 

m 

c c 

c 

c 

c 

C 

c 

c 


w-« H 


(V! 

(Vi OJ 

OI OI 

cc 

o 

C 

C CJ 

o o 

♦ 4 

♦ 

♦ 

♦ 

4 4 

UJ uJ 

UI 

ul 

Ui Ui 

UI UI 

in tv 

— 

o. 

>o re 

in 4 


(VI 

4 

so rr 

r< .r 

rr 4 

<v 

CT 

fv (T 

4 c 

<r o 

Oi 

>o 

0 Oi 

<r ao 

<r»4 

w* 

vO 

<r v0 

— fv 

mn 

o> 

OJ 

>o *-• 

nq in 

re — 

rv 

c 

«o rv 

>c cr 

>c cr 

<r 

o 

OJ >0 

— >0 

4 X 

rv 

*- 


o o 

• • 

• 

• 

• » 

• • 

oo 

o 

o 

o c 

o o 


oo 

o 


c o 

c o 

Co 

o 

O 

o r j 

o o 

♦ 4 

4 

♦ 

♦ 

4 

♦ ♦ 

uiUI 

UI 

UI 

UI Ul 

ul UJ 

rv >o 

o. 

o 

fv n 

G X 

som 

Oj 

ac 

in 

re 

<i ao 

nC 4 

X 

<r 

fv fv 

fv O' 

O' *a 

*C 

Oj 

4 

re 

tr rv 

ere 

C 

X 

4 

rv 

4 rr 

4 4 

<3 

«-• 

ao 

rv 

4 4 

-*o 

*-» 

ao 

30 

O 

OI >o 

in g 

4 

Oi 

O 

X 

4 00 

*- rr 

4 

IT 

>£ 


rv rv 

• * 

» 

• 

• 

• 

• • 


U 

O 

O 

o 

O CJ 


c r: 

o 

c 

o 

o 

o o 

*-* g 

o 

*_> 

o 

g 

u o 

♦ ♦ 

♦ 

♦ 

4 

4 

4 4 

LxJ W 

Ul 

ui 

Ul 

UJ 

Ul Ul 

rv rr 

4 

Oj 

re 

41 

fv a\ 

X3 — 


rr 

4 

X 

) rv 

O' O' 

X 

cr 

re 

X 

re x 

0>G\ 

C\ 

cr. 

On 

oi 

oi a. 

T 

4 

;r 

rv 

rv 

er x 

G JN 

<r 

(V 

rv 

X 

4 O' 

rr — 

*- 

C 



cr (r 

(X c 

c 

*— 

ex 

re 

in x 

(VJOJ 

Oi 

OJ 

OJ 

OJ 

Oj oj 

• • 

• 

• 

• 

• 

• • 

c c 

c 

c: 

c 

c 

c c 


n 4 

4 

4 

4 

m 

r. 

c e 

G 

rr 

o 


C 'J 

4 4 

4 

4 

4 

4 

4 4 

UJ u 

Ul 

LU 

w 

Ul 

U jJ 

x evj 

4 

a0 

4 

m 

X IT 

rv Oj 

X 

fv 

*-♦ 


ry 4 

n cr 

4 

X 

a* 


O' C 

O. 4 


O' 



c *-< 

«m 

Oj 

4 


4 

f*) fv 

4) N 

rv 

CVI 

n 

O 

a\ n 

4 4) 

♦ 

OJ 

o 

OJ 

O 4 

OjjO 

*-* 

n 


in 

n *-• 

OlO> 

4 

Oi 

«-« 

rv 

4 ™ 

• • 

• 

• 

• 

• 

• • 

C G 

C 

o 

G 

c 

o c 



4 >0 

00 o> 


a 0 >0 4 OI 

o^tum 


O 0 

4 * 


D-28 



8*»* (!•«!( -'ll 99r .^ *(1 n»17B101395E*r<3 C • 3 1 635 U73E +00 

9 fc» n*?ibl 1 -M, J !"!• 16 . 8 <**f< 4 rtCE*r ;3 n • ? 1 83 «P‘ib 7 E + 00 

108* 0.2bt6o*'v:j9 i >r - n*16?C6S>8*7f 403 0177<*40*9f,k +C0 


c c c c e 

O O O O CJ 


Ui oi UJ Uj oJ 

dKhiu 
C\i>Cirr' r * 
ac rv cv cr cr 

0 C (V. (T « N 

a ~ cr c * 
ownrvcvi* 

^ in r« in *n 

vc o oc u 
— tv 



C C O C G 


<r r <*•, w it 
g; c o u o 
♦ ♦ ♦ ♦ ♦ 
Ll UJ u; it aJ 
IT ^ r* a 'J 
'T <*+ (\* «-* lf> 

x ac 4- rr \p 
o rtj uc c n 

at * (V — rv 

^ rs oc 
r* <*i <r tr tv. 
x * 4 4 - 



c c c c c 



c c. r: c c 

♦ ♦ ♦ + ♦ 

aJ 4 ' l/J> 'k ij 

foj) ^ *» (T> 
^'0 4 r r 
rv *jc r x r* 
O' 0 ‘ r * a> 
N JC 
— /I 0 A 4 t 4 > 
<f\ 4 If) 4 

c* sC ^ rn ui 
.v n 4 nC t) 
• * * • « 
c o c o c 




D-29 



84. 0«JC4<'61O9F4Cfi -n*877964n?Ct+C0 fi.4392*3KbfiE 4fi? -fi.1176r4i79E.4-c2 

9fc. 0*4224 7fi927e *fifi fi . ?379f>3 1 26E 4tfi C . 1 1 42fib 1 1 ?E 4fi? 0 . 1 1 4*- C3369E 4fi 1 

106. fi . b3 1 £C 70 45E *fin fi» 194M Ch96E4Cl C . S474«i3n??E 4fil 0.«fr5948*58?E4Cl 


nj <u ru (U 
C- O O Q O 
♦ 4 ♦ ♦ ♦ 

uJ LU L*J LkJ (L*J 

'C a ^ n r 
to :j a n 

N 4 <t r 4 

n m rvj — » m 

rr ff\ n -* <v 

IT t O flj 4 

(*IIT O' « AC 
— ca- 4 a) 



c o o c c 


a o u o o 

♦ ♦ ♦ ♦ ♦ 

JjJ lU uJ LjJ UJ 

m cr: iv. sC 
O ip at) O lO 

4 cr rv 4 c 
h n n at iC 
o: a rr x oc 
rv fv ^ fv ip 

m in ru ao o 

O' Oil h o> 

(Ttrrvcv*- 

• • • • • 

o c c u c 



^ *• aj (v 

O t- u o o 

+ ♦ ♦ ♦ + 

'.ii ui Ui u.‘ u 

m c\j oc iu ru 

fv /) v u N r> 
-vC4NO 
CU 4 0 s o 
x in vC 4© 
C -X 

x o. x rr cr 
IT X rv 0 j it 
t n oc 4 


c c c c c 



a o c c a 

c o a c o 

+ + + + + 

u. il i) U aJ 
fT f\j 'T Oj CT^ 
- t C fV 4 

rv x (\* ic n 

P* a", in o >£< 

X >o ru * 
iO P4. .u o ru 
n 4 ru rv so 
o> o 4 oi 

n in n n 4 



cocoa 


4 .0 o >n 
HOjm 4 4 



D-30 



CitlG2CiJ3Cnbt4.cn *n»l t )3?91?e3E*C. 1 - C. • 4386*39E'4E ♦ <;? 

n»9lH36bfe57t400 f!» 147554 1 1 1E + 01 C. • A&599 1 1 9?E *C. 1 


fV fUV. (V (V fv. 

ooca(;c 

♦ > ♦ ♦ ♦ ♦ 

1 1 . jxl LlI ui L *■ ■ Laj 

a; tv a ^ u x 

^ m u ui fv ai 

>c «- * r* rv 

^ -fi -» * * <. 

cm ^ x a 
>o >o -t n ^ 

^ vO fv rvj 

oc r\ <r* m o in 

«-aa <r ^ 4 

• ••••• 

c c c c c a 



c o u c o u 

♦ ♦♦♦♦♦ 

UJ UJ -xj Lu iaJ uJ 

0> PC 0> C\j (J\ 

M *0^0 O 
«~i rc r* p*. *~ 
ir> *• o * ^ 
no(V 4 aa c 
(v ■* o a a •* 

(X. * ir c ir <r 

-» iP U'i * JO 

<\u cvj ru t\j (\j 



c c c c cz c 


c o c o c c 
a o a o c c 

♦ ♦♦♦♦♦ 

-J .jJ -»J jJ ui 

CVi fV O <T> r-t IT5 

h ^ ^ n 

♦ * X ^ 30 (V 

(T U i OC n 

^ jo rx *** no >0 
H IX (1) ,VJ ♦ 

>0 ^ (V O' M x 

rvj ^ (v rx * *-» 

X) N N vO X >0 

u c c c o a 




a 

T> 

r- 

Ji 


D-31 



